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INTRODUCTION 


Numerous  investigations  have  been  carried  out  during  the  past 
decade  into  the  application  of  synthetic  aperture  radar  (SAR)  technol¬ 
ogy  to  various  oceanographic  measurement  and  detection  problems. 
These  investigations  utilized  the  Seasat  SAR  and  several  aircraft- 
mounted  SAR  systems,  which  differ  in  many  details  but  are  all 
representative  of  the  same  general  type  of  system,  i.e.,  the 
conventional  strip-map  SAR  which  is  described  in  Section  2  of  this 
report.  The  results  of  these  investigations  indicate  that  strip-map 
SAR  systems  can  provide  useful  information  on  a  variety  of  phenomena 
including  surface  gravity  waves,  internal  waves,  bottom  topographic 
features,  ocean  currents,  fronts,  winds,  surface  films,  and  sea  ice. 
However,  these  systems  are  also  subject  to  certain  limitations  and 
restrictions.  These  include  operation  at  a  fixed  look  direction  and  a 
single  transmitted  wavelength.  In  an  attempt  to  relieve  these 
limitations,  several  modifications  of  the  conventional  SAR  system 
design  have  been  proposed  in  recent  years,  some  of  which  have  been 
implemented  and  some  of  which  are  still  conjectural.  These  advanced 
SAR  concepts  are  described,  and  their  possible  advantages  for  ocean 
remote  sensing  are  discussed  in  Section  3  of  this  report. 

In  order  to  evaluate  the  utility  of  SAR  for  measuring  surface 
gravity  wave  spectra,  as  well  as  for  other  ocean  applications,  a  nu¬ 
merical  image  simulation  model  has  been  developed  and  tested  using 
SIR-B  data.  This  model  is  described  in  Appendix  A  of  this  report. 
Additionally,  efforts  have  been  made  In  collaboration  with  F.  Monaldo 
of  the  Applied  Physics  Laboratory  to  develop  procedures  for  estimating 
wave  height  and  slope  spectra  from  SAR  Images.  A  published  report  on 
the  results  of  this  collaboration  are  reproduced  in  Appendix  B  of  this 
document. 
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CONVENTIONAL  SAR  SYSTEMS 


The  conventional  type  of  SAR  system,  which  Includes  Seasat,  SIR-A, 
and  SIR-B  as  well  as  several  different  aircraft  SAR  systems,  consists 
of  a  transmitter  which  generates  a  sequence  of  microwave  pulses,  an 
antenna  which  Is  maintained  In  a  fixed  orientation  relative  to  the 
flight  direction,  a  receiver  which  measures  and  records  the  amplitude 
and  phase  of  the  returned  signals,  and  a  processing  system  which  re¬ 
constructs  the  recorded  signals  Into  a  map  or  image  of  the  scene. 
This  type  of  system  is  commonly  referred  to  as  strip-map  SAR.  A  dia¬ 
gram  of  the  data  collection  geometry  for  this  type  of  system  Is  shown 
in  Figure  1. 

The  transmitted  pulses  are  usually  In  the  form  of  linear  FM  sig¬ 
nals  which  are  converted  Into  equivalent  short  pulses  by  convolving 
the  returned  signal  with  a  reference  signal  having  a  form  similar  to 
that  of  the  transmitted  pulse.  In  this  way,  a  high  spatial  resolution 
In  the  cross-track  or  range  direction  Is  obtained.  The  range  resolu¬ 
tion  in  this  type  of  system  Is  given  by 


where  c  Is  the  speed  of  light  and  Br  is  the  bandwidth  of  the  transmit¬ 
ted  signal . 

In  the  along-track  or  azimuth  direction,  a  similar  method  is  used 
except  that  the  frequency  modulation  Is  obtained  from  the  Doppler 
shift  of  the  reflected  signals  due  to  the  platform  motion.  The  azi¬ 
muth  resolution  Is  given  by 


=  Ba 
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where  V  is  the  platform  velocity  and  Ba  Is  the  azimuth  bandwidth, 
which  Is  related  to  the  antenna  beamwidth  p,  the  platform  velocity  V, 
and  the  radar  wavelength  X  by  the  equation 

Ba=M.  (3) 

Since  the  antenna  beamwidth  is  given  by  the  ratio  of  the  radar  wave¬ 
length  to  the  antenna  length,  the  azimuth  resolution  is  nominally 
equal  to  half  the  antenna  length  [1]. 

In  actuality,  Eq.  (2)  represents  a  limiting  value  for  the  resolu¬ 
tion,  which  is  not  attained  if  there  are  phase  errors  due  to  scatterer 
motions  or  variations  In  the  platform  velocity.  Nevertheless,  the 
azimuth  resolution  of  a  SAR  Is  improved  (within  limits)  by  broadening 
the  antenna  beam  rather  than  by  narrowing  it  as  In  a  real  aperture 
radar  system.  Also,  the  azimuth  resolution  is  independent  of  range 
Instead  of  being  proportional  to  range  as  In  a  real  aperture  system. 
Thus,  the  synthetic  aperture  technique  allows  fine  resolution  imagery 
to  be  obtained  from  high  altitudes.  For  example,  to  obtain  the  same 
resolution  as  the  Seasat  SAR  with  a  real  aperture  radar  at  the  same 
altitude  and  wavelength  would  require  an  antenna  more  than  15  km  long. 

The  Increased  resolution  obtainable  by  synthetic  aperture  tech¬ 
niques  is  obtained  at  the  cost  of  a  larger  amount  of  processing  than 
is  required  for  real  aperture  radars.  Usually,  this  processing  Is 
carried  out  as  a  separate  stage,  after  the  data  collection  is  com¬ 
pleted.  However,  developments  in  digital  processing  technology  have 
made  real-time  on-board  processing  feasible,  and  SAR  systems  utilizing 
such  processors  now  exist  [2],  Research  is  continuing  in  the  develop¬ 
ment  of  more  efficient  processors  as  well  as  algorithms  for  the  ex¬ 
traction  of  specialized  information  from  SAR  data  [3]. 

Because  the  azimuth  resolution  in  a  SAR  is  obtained  by  utilizing 
the  Doppler  shift  of  the  returned  signal  as  the  object  passes  through 
the  antenna  beam,  the  imaging  process  is  Influenced  by  the  motion  of 
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the  objects  imaged.  An  object  having  a  constant  velocity  Vr  in  the 
line  of  sight  direction  is  displaced  in  the  image  by  an  amount 

Ax  =  |  Vr  (4) 

in  the  along-track  direction,  where  R  is  the  radar  range  and  V  is  the 
platform  velocity.  If  the  object  changes  its  radial  velocity  ( i . e . , 
accelerates)  during  the  time  it  is  within  the  antenna  beam,  its  image 
is  smeared  in  the  along-track  direction  by  an  amount  proportional  to 
the  change  in  its  radial  velocity  during  this  time. 

In  the  case  of  the  ocean  surface,  the  Image  displacements  caused 
by  the  motion  of  the  surface  can  cause  variations  in  the  image  inten¬ 
sity  which  may  either  enhance  or  degrade  the  visibility  of  surface 
features  such  as  gravity  waves  or  current  boundaries.  The  orbital 
motions  associated  with  gravity  waves  can  cause  these  waves  to  be 
imaged  by  a  SAR  because  of  the  differential  displacement  (or  velocity 
bunching)  effect,  even  when  the  waves  would  not  be  imaged  by  a  real 
aperture  radar  system  [4].  The  effects  of  scatterer  accelerations  and 
random  motions,  on  the  other  hand,  cause  a  degradation  in  the  azimuth 
resolution  which  may  smear  out  a  pattern  which  would  otherwise  be 
resolvable  [5]. 

The  following  paragraphs  summarize  the  current  status  of  investi¬ 
gations  into  the  use  of  SAR  for  various  oceanographic  applications. 
This  summary  is  not  comprehensive,  but  is  intended  to  give  an  overview 
of  this  field  of  research  as  well  as  some  recent  results. 

2.1  SURFACE  GRAVITY  WAVES 

The  oceanographic  application  of  SAR  which  has  received  perhaps 
the  most  attention  to  date  is  the  measurement  of  surface  gravity  wave¬ 
lengths  and  directions.  A  number  of  field  experiments  dedicated  to 
this  problem  were  conducted  during  the  1970's,  including  the  1975 
Marineland  experiment  [6],  several  experiments  held  during  the  Seasat 
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mission  [7,  8],  and  the  1979  MARSEN  experiment  [9].  More  recent  ex¬ 
periments  such  as  those  conducted  during  the  SIR-B  mission  [10]  have 
also  addressed  this  problem  using  a  wider  range  of  wave  conditions, 
different  sensor  parameters,  and  more  accurate  i_n  situ  measurements. 

The  experiments  conducted  to  date  have  demonstrated  that  under 
certain  environmental  conditions  and  SAR  system  configurations,  the 
dominant  wavelength  and  direction  can  be  accurately  estimated  from  SAR 
data.  The  precise  range  of  conditions  under  which  these  quantities 
are  measurable  is  a  complicated  function  of  the  SAR  system  as  well  as 
environmental  parameters,  and  has  not  yet  been  completely  defined. 
However,  in  general  waves  which  are  traveling  in  the  range  direction 
(see  Figure  1)  are  more  easily  detected  and  more  accurately  measured 
than  those  traveling  in  the  along-track  direction  [11-12]. 

An  example  is  provided  by  the  Mari nel and  SAR  data  set  which  was 
collected  with  eight  different  look  directions  on  successive  over¬ 
flights.  Figure  2  shows  the  X-  and  L-band  images  from  three  of  these 
passes,  for  which  the  look  direction  was  approximately  cross-wave 
(upper  images),  down-wave  (middle)  and  up-wave  (lower  images).  Visual 
examination  of  these  images  indicates  that  range-traveling  waves  are 
more  easily  discernible  than  azimuth-traveling  waves.  A  quantitative 
measure  of  the  wave  contrast,  or  modulation  depth  (defined  in  [11]) 
also  indicates  this  dependence,  as  shown  in  Figure  3.  Furthermore, 
estimates  of  the  dominant  wavelength  and  direction  obtained  from  each 
pass  are  more  accurate  for  the  range  traveling  cases  than  the  azimuth¬ 
traveling  cases,  as  shown  in  Figures  4  and  5  (taken  from  Ref.  [12]). 

In  addition  to  the  experimental  studies  mentioned  above,  a  great 
deal  of  theoretical  work  has  been  carried  out  regarding  the  mechanism 
by  which  SAR  images  gravity  waves,  in  order  to  predict  the  range  of 
detectability  of  these  waves  and  the  accuracy  of  their  measurement  by 
SAR.  Numerous  papers  have  been  written  on  this  subject,  and  an  aura 
of  controversy  has  surrounded  the  problem.  A  recent  review  and  re¬ 
formulation  of  this  problem  by  Hasselmann,  et  al .  [13]  reflects  the 
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Figure  3.  Modulation  Depth  Versus  SAR  Look  Direction 
for  Marineland  Data  Set  [11] 


Figure  4.  Apparent  (SAR-Derived)  Wavelength  Versus  SAR 
Look  Direction  for  Marineland  Data  Set  [12] 
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Figure  5. 


Apparent  (SAR-Derived)  Wave  Direction  Versus 
SAR  Look  Direction  for  Marine! and  Data  Set 
[12] 
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opinions  of  many  investigators  although  there  is  still  disagreement 
regarding  fundamental  issues  as  well  as  the  details  of  the  imaging 
process  [14,  15]. 

Because  of  the  complexity  of  the  SAR  imaging  process  as  well  as 
the  surface  description,  numerical  simulation  models  have  been  found 
to  be  useful  for  understanding  and  predicting  the  SAR  response  to 
gravity  waves  [16,  17].  In  the  simulation  model  implemented  at  ERIM, 
the  surface  is  assumed  to  be  described  by  a  wave  height  spectrum  with 
statistically  independent  spectral  components.  A  particular  realiza¬ 
tion  of  the  surface  consistent  with  this  description  is  computed  by 
assigning  random  phases  to  each  spectral  component.  The  radar  cross 
section  of  the  surface  is  computed  for  each  resolution  cell  using  a 
modulation  transfer  function  which  Incorporates  both  hydrodynamic  and 
tilt  effects  [4],  and  the  surface  motion  effects  are  included  by  com¬ 
puting  the  relevant  motion  parameters  for  each  resolution  cell  [17]. 

Example  results  of  this  simulation  procedure  are  shown  in 
Figure  6.  These  images  represent  the  predicted  response  of  the  Seasat 
and  SIR-B  systems  to  the  wave  height  spectrum  shown  in  Figure  7,  which 
has  a  significant  wave  height  of  3.6  meters  and  a  dominant  wavelength 
of  250  meters.  The  direction  of  the  waves  for  the  three  cases  shown 
in  Figure  6  are  0°,  45°  and  90°  with  respect  to  the  radar  look  direc¬ 
tion.  The  results  for  Seasat  confirm  the  experimental  findings  that 
gravity  waves  are  more  sensitively  detected  and  more  accurately  repre¬ 
sented  on  SAR  images  for  which  the  look  direction  is  nearly  parallel 
to  the  wave  propagation  direction,  than  on  those  for  which  the  look 
direction  is  at  a  large  angle  to  the  wave  propagation  direction.  In 
particular,  the  apparent  wave  direction  on  the  Seasat  image  for  the 
45°  case  is  shown  in  Figure  8  to  deviate  toward  the  range  direction 
from  the  actual  wave  direction,  by  approx imate i /  8°  ‘or  this  case. 

The  results  for  SIR-B  show  that  the  SAR  oaqes  and  spectra 
(Figure  9)  much  more  accurately  represent  the  actual  wave  conditions 
than  the  Seasat  images.  The  difference  is  due  to  the  fact  that  the 
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Figure  6. 


Seasat  and  SIR-B  Simulations  for  Wave 
Conditions  Described  By  Spectrum  Shown 
in  Figure  7,  for  Three  Look  Directions 
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Figure  7.  Wave  Height  Spectrum  Used  in  Simulations,  Expressed 
as  a  Function  of  (a)  Wave  Number,  and  (b)  Wave 
Direction 
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Figure  8. 
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Seasat  Image  Spectrum  for  45°  Look  Direction, 
Reduced  to  1-D  Spectra  in  Wave  Number  and 
Angle 


R/V  ratio  [see  Eq.  (4)]  is  equal  to  approximately  120  sec  for  Seasat 
and  35  sec  for  SIR-B.  Thus,  the  degradation  in  image  quality  due  to 
scatterer  motions  is  much  less  severe  for  SIR-B  than  Seasat.  These 
results  have  been  confirmed  by  analyzing  the  data  collected  during  the 
SIR-B  mi ssion  [10,  17] . 

The  implication  of  these  results  Is  that  the  look-direction  limi¬ 
tations  inherent  in  strip-map  SAR  systems  can  be  at  least  partially 
obviated  by  the  use  of  platforms  with  smaller  R/V  ratios.  Another 
possible  solution  to  this  problem  is  the  use  of  angle  diversity  or 
spotlight  SAR  systems  which  allow  the  look  direction  to  be  varied  from 
the  fixed  (broadside)  direction  used  in  strip-map  SARs.  This  concept 
is  further  explored  in  Section  3  of  this  report. 

2.2  INTERNAL  WAVES 

Oceanic  Internal  waves  have  been  observed  on  both  aircraft  and 
Seasat  SAR  images  [18-20]  and  have  recently  been  the  subject  of  sev¬ 
eral  dedicated  field  experiments  [21,  22].  The  principal  mechanism 
responsible  for  the  imaging  of  internal  waves  by  SAR  is  thought  to  be 
the  interaction  of  short  surface  waves  with  the  Internal  wave-induced 
surface  currents.  However,  the  damping  effects  of  surface  films, 
which  are  alternately  compacted  and  dissipated  by  the  internal  wave 
surface  currents,  may  also  be  a  factor.  Model  calculations  based  on 
the  former  mechanism  appear  to  slightly  underpredict  the  observed 
features  at  L-band  and  to  grossly  underpredict  the  patterns  observed 
on  X-band  imagery. 

Although  further  research  appears  to  be  needed  to  quantitatively 
understand  the  contrast  of  the  observed  internal  wave  features  (par¬ 
ticularly  at  X-band),  the  existence  and  form  of  these  features  provide 
useful  oceanographic  information  about  the  generation  and  propagation 
of  internal  waves.  The  presence  of  these  features  also  yields  infor¬ 
mation  on  stratification  conditions  and  possibly  on  bottom  topographic 
features  which  cause  the  internal  waves  to  be  formed. 


The  primary  limitation  of  strip-map  SAR  systems  for  observing 
internal  waves  is  their  sensitivity  to  look  direction.  The  scattering 
of  microwave  radiation  from  the  surface  is  controlled  primarily  by  the 
amplitude,  or  spectral  density,  of  waves  traveling  in  the  range  direc¬ 
tion  and  having  wavelengths  which  satisfy  the  Bragg  criterion,  i.e., 

XB  =  Xq/(2  sin  9)  (5) 

where  X0  is  the  radar  wavelength  and  9  is  the  incidence  angle.  These 
waves  interact  with  the  surface  currents  induced  by  internal  waves, 
and  the  resulting  variations  in  their  amplitude  or  spectral  density 
cause  variations  in  the  surface  reflectivity,  which  are  detected  by 
the  SAR.  The  strength  of  the  wave-current  interaction  is  highly  de¬ 
pendent  on  the  relative  angle  between  the  current  direction  and  the 
Bragg  wave  propagation  direction  (or  radar  look  direction).  In  par¬ 
ticular,  internal  waves  propagating  perpendicularly  to  the  line  of 
sight  direction  cause  smaller  changes  in  the  Bragg  wave  amplitude,  and 
therefore  are  less  detectable,  than  those  traveling  in  the  line-of- 
sight  direction. 

In  the  conventional  strip-map  SAR  configuration,  the  look  direc¬ 
tion  is  always  normal  to  the  ground  track.  As  a  result,  internal 
waves  traveling  parallel  to  the  ground  track  are  less  detectable  than 
those  traveling  perpendicularly  to  the  ground  track.  Angle  diversity 
techniques  remove  this  limitation  by  allowing  the  look  direction  to  be 
varied  (in  spotlight  mode)  or  by  using  two  or  more  fixed  look  direc¬ 
tions  (in  squint  mode),  for  example  at  angles  of  *45  to  the  normal 
broadside  viewing  direction. 

Additionally,  since  the  strength  of  the  wave-current  interaction 
depends  on  the  wavelength  of  the  surface  waves,  the  detectability  of 
Internal  waves  may  be  optimized  by  using  the  frequency  diversity  tech¬ 
niques  discussed  in  Section  3.2.  These  may  Involve  either  changing 
the  frequency  of  the  transmitted  signal,  and  thus  varying  the  Bragg 


wavelength,  or  using  two  closely-spaced  frequencies.  The  latter  tech¬ 
nique  effectively  causes  the  system  to  respond  to  surface  waves  with 
wavelengths  which  are  determined  by  the  difference  frequency  [see 
Eq.  (9)]. 

2.3  BOTTOM  TOPOGRAPHY 

Several  types  of  features  on  SAR  images  can  be  used  to  make  infer¬ 
ences  regarding  bottom  topography  [23].  These  include  internal  wave 
patterns  (as  discussed  in  the  previous  section),  surface  gravity  waves 
with  wavelengths  sufficiently  long  to  interact  with  the  bottom,  and 
patterns  associated  with  local  variations  in  tidal  currents  due  to  the 
interaction  of  these  currents  with  the  bottom  topography. 

The  latter  type  of  bottom-related  features  show  a  striking  amount 
of  detail.  Well-known  examples  include  Seasat  images  of  the  English 
Channel  and  the  Nantucket  Shoals.  Several  modeling  studies  using  an 
approach  similar  to  that  used  in  the  modeling  of  Internal  wave  signa¬ 
tures  have  been  carried  out  [24-26].  These  models  yield  generally 
good  agreement  with  observed  L-band  SAR  intensity  modulations  except 
for  certain  geometries  which  produce  much  larger  modulations  than 
predicted  by  the  models. 

Because  of  the  relative  success  in  modeling  these  features,  there 
is  some  hope  of  developing  an  inversion  technique  to  calculate  water 
depths  from  the  observed  SAR  signals.  However,  since  these  signatures 
are  highly  dependent  on  the  wind  conditions  as  well  as  the  currents,  a 
knowledge  of  the  wind  and  current  speeds  and  directions  would  be 
needed  to  carry  out  such  an  inversion. 

The  SAR  signatures  associated  with  bottom  topography  are  also 
look-angle  dependent  and  therefore  could  be  optimized  by  the  use  of 
angle  diversity  techniques.  The  optimum  look  direction  is  in  general 
parallel  to  the  direction  in  which  the  gradient  of  the  water  depth  is 
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a  maximum.  Effects  of  scatterer  motions  are  not  critical  to  the  for¬ 
mation  of  these  signatures,  although  the  resolution  degradation  due  to 
random  surface  motions  may  contribute  to  the  look-angle  dependence. 

2.4  OCEAN  CURRENTS 

Current  boundaries  are  frequently  visible  on  SAR  images,  possibly 
because  of  variations  in  surface  roughness  as  well  as  surface  velocity 
across  the  boundary.  Variations  in  surface  roughness  are  manifested 
as  a  change  in  the  image  intensity  (tone)  across  the  boundary,  while 
changes  in  the  surface  velocity  manifest  themselves  through  the  Dop¬ 
pler  shift  effect  [27],  These  image  patterns  may  be  used  to  chart  the 
position  of  large-scale  current  features  such  as  the  Gulf  Stream. 

A  quantitative  measurement  of  the  surface  current  speed  may  also 
be  possible  by  performing  a  spectral  analysis  of  the  recorded  SAR 
signals  or  the  complex  SAR  image  [28,  29],  The  Doppler  spectrum  of 
these  signals  is  determined  by  the  antenna  pattern,  the  platform  ve¬ 
locity  and  the  radar  wavelength.  For  stationary  surfaces,  the  maximum 
spectral  density  occurs  at  zero  frequency,  assuming  the  antenna  is 
aligned  perpendicularly  to  the  flight  path.  If  the  entire  surface 
within  the  field  of  view  of  the  antenna  is  moving  with  a  radial  veloc¬ 
ity  Vr,  the  Doppler  spectrum  is  shifted  by  an  amount 


where  X  is  the  radar  wavelength.  If  this  shift  is  an  appreciable 
fraction  of  the  Doppler  bandwidth  Ba,  given  in  Eq.  (3)  above,  and  if 
the  antenna  pointing  angle  remains  stable,  the  shift  can  be  measured 
and  used  to  infer  the  radial  component  of  the  surface  velocity.  Mea¬ 
surements  of  surface  currents  in  the  range  of  1-2  m/s  have  been  mea¬ 
sured  using  aircraft  X-band  SAR  data  in  this  way  [30,  31]. 

Attempts  to  apply  this  method  to  Seasat  data  have  not  been  suc¬ 
cessful  because  the  Doppler  bandwidth  due  to  the  platform  motion  is 
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too  large  compared  to  the  Doppler  shifts  induced  by  surface  currents. 
The  Doppler  bandwidth  can  be  reduced  by  narrowing  the  antenna  beam- 
width,  as  shown  by  Eq.  (3),  but  this  implies  a  degradation  in  azimuth 
resolution,  for  a  given  radar  wavelength.  Comparison  of  Eqs.  (2),  (3) 
and  (6)  shows  that  for  Af/Ba  to  be  on  the  order  of  0.1,  while  main¬ 
taining  a  resolution  />a,  the  radar  wavelength  must  be  on  the  order  of 

X  ~  20 />aVr/V  (7) 

in  order  to  measure  a  radial  velocity  Vr.  This  analysis  indicates 
that  an  X-band  SAR  having  the  same  resolution  as  Seasat  and  operating 
at  the  same  altitude  should  be  marginally  able  to  detect  a  radial 
velocity  of  1  m/s.  Measurements  of  smaller  currents  would  require  a 
narrower  antenna  beam  and  consequently  a  coarser  azimuth  resolution. 

An  alternative  method  of  measuring  surface  currents  using 
multiple-antenna  SAR  systems  is  suggested  in  Section  3.3.  This  method 
could  also  be  combined  with  angle  diversity  methods  to  measure  both 
components  of  the  current. 

2.5  WINDS 

Since  the  radar  cross  section  of  the  ocean  surface  is  dependent  on 
the  wind  speed,  it  is  possible  to  determine  the  wind  speed  and  direc¬ 
tion  by  a  series  of  radar  cross  section  measurements  taken  with  dif¬ 
ferent  look  directions  [32].  The  same  principle  allows  an  estimate  of 
wind  speed  to  be  made  from  the  average  image  intensity  within  a  SAR 
image  [33]  although  the  wind  direction  cannot  be  determined  by  this 
method  since  only  one  look  direction  is  usually  employed.  Again, 
angle  diversity  methods  offer  the  possibility  of  overcoming  this 
1  imitation. 

A  second  method  of  obtaining  wind  information  from  SAR  images  is 
to  examine  the  imagery  for  the  existence  of  wind  rows,  or  roll  vorti¬ 
ces  [34].  These  features  can  be  quantified  by  doing  a  spectral  analy¬ 
sis  of  the  image,  and  the  results  used  to  estimate  the  wind  speed  and 
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direction.  The  spacing  of  the  rows  is  related  to  the  wind  speed  and 
their  orientation  may  be  used  to  infer  the  wind  direction  [35].  These 
features  also  appear  to  be  look  angle  dependent  and  would  benefit  from 
the  use  of  multiple  look  directions. 
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NON-CONVENTIONAL  SAR  SYSTEMS 


Several  variations  on  the  conventional  SAR  system  design  have  been 
either  proposed  or  constructed  in  recent  years  [36].  These  non- 
conventional  systems  are  described,  and  their  potential  applicability 
to  oceanographic  measurement  problems  are  discussed  in  the  following 
sections. 

3.1  ANGLE  DIVERSITY  TECHNIQUES 

The  along-track  resolution  of  a  SAR  is  determined  by  the  Doppler 
bandwidth  of  the  signals,  which  is  in  turn  determined  by  the  range  of 
angles  over  which  a  given  object  is  illuminated.  In  a  conventional  or 
strip-map  SAR,  this  range  of  angles  is  equal  to  the  antenna  beamwidth, 
since  the  antenna  is  held  in  a  fixed  orientation  perpendicular  to  the 
flight  direction.  Several  alternative  modes  of  operation,  referred  to 
collectively  as  angle  diversity  techniques,  have  been  devised  in  which 
the  antenna  is  not  held  in  this  fixed  orientation. 

In  one  of  these  modes,  referred  to  as  spotlight  mode,  the  antenna 
is  continuously  rotated  so  as  to  illuminate  a  given  area  on  the  sur¬ 
face.  Other  variations  include  squint  mode,  in  which  the  antenna  is 
held  in  a  fixed  but  arbitrary  orientation,  and  swath  mode,  in  which 
the  antenna  orientation  is  programmed  to  map  out  an  arbitrary  path  on 
the  surface.  A  diagram  of  the  data  collection  geometry  for  these 
modes  is  shown  in  Figure  10. 

In  spotlight  mode,  the  along-track  resolution  is  controlled  by  the 
angle  hd  over  which  the  antenna  is  rotated,  and  the  scene  size  is 
determined  by  the  antenna  beamwidth.  In  order  to  achieve  this  resolu¬ 
tion,  however,  a  modified  method  of  recording  and  processing  the  data 
is  required.  The  polar  format  technique  developed  by  Walker  [37] 
makes  the  spotlight  mode  of  operation  feasible.  This  technique  has 
been  implemented  optically  by  recording  the  time  histories  for  each 


returned  pulse  along  a  series  of  radial  lines  on  a  circular  film, 
rather  than  along  a  series  of  parallel  lines  on  a  strip  of  film  as  is 
done  in  the  conventional  SAR  case.  In  the  digital  implementation,  the 
same  result  is  obtained  by  a  resampling  procedure. 

During  processing,  the  entire  data  set  can  be  processed  coherently 
in  order  to  maximize  the  resolution,  or  the  data  can  be  processed  in 
segments  at  a  coarser  resolution  in  order  to  obtain  a  number  of  inde¬ 
pendent  looks  or  images  of  the  scenes.  For  a  segment  corresponding  to 
an  angular  Interval  A^i ,  the  azimuth  resolution  is  given  by 

pi  =  ^ 

where  X  is  the  radar  wavelength.  If  N  such  angle  Intervals  are  pro¬ 
cessed  and  Incoherently  summed,  the  signal  variance  due  to  speckle  on 
a  distributed  scene  is  reduced  by  a  factor  of  1/N  over  that  for  a 
fully  coherent,  or  single-look  image. 

Other  methods  of  processing,  or  information  extraction,  may  be 
beneficial  if  the  scene  Imaged  is  either  time-dependent  or  aspect- 
angle  dependent.  These  methods  would  make  use  of  the  fact  that  data 
can  be  collected  by  angle  diversity  systems  over  an  extended  time 
interval  (a  few  seconds  to  a  few  tens  of  seconds)  as  well  as  an  ex¬ 
tended  angular  interval.  For  example,  gravity  wave  images  collected 
over  an  extended  period  of  time  could  be  used  to  obtain  independent 
spectral  estimates  which  can  be  averaged  together  to  Increase  the 
number  of  degrees  of  freedom.  Alternatively,  images  formed  from  each 
sub-interval  could  be  compared  with  each  other  to  look  for  differ¬ 
ences,  rather  than  summed  together  to  reduce  speckle. 

Differences  in  time  between  sub-intervals  might  be  used  to  quan¬ 
tify  ice  movements  or  surface  currents,  or  to  determine  the  direction 
of  propagation  of  gravity  waves  and  thus  resolve  the  180°  ambiguity 
Inherent  in  conventional  SAR  measurements  of  wave  direction.  Differ¬ 
ences  in  look  direction  between  sub-intervals  could  be  used  to 


optimize  the  detectability  of  surface  patterns  which  have  a  strong 
angle  dependence,  or  to  classify  surface  features  on  the  basis  of  this 
dependence.  Similar  advantages  may  apply  to  squint  mode  systems, 
particularly  those  employing  multiple  antennas,  as  described  in 
Section  3.3  of  this  report. 

3.2  FREQUENCY  DIVERSITY  TECHNIQUES 

Another  category  of  advanced  SAR  system  concepts  concerns  the 
utilization  of  multiple  radar  frequencies  or  frequency  intervals.  One 
type  of  multi -frequency  SAR  Is  that  which  collects  data  simultaneously 
but  not  coherently  at  two  or  more  widely  separated  wave-lengths.  An 
example  of  an  existing  SAR  system  of  this  type  is  the  ERIM/CCRS  X-L-C 
band  SAR,  which  collects  data  simultaneously  at  either  X  and  L  band  or 
X  and  C  band. 

The  advantage  of  this  type  of  system  is  that  it  allows  information 
about  ocean  processes  and  scattering  phenomena  to  be  inferred  from  the 
wavelength  dependence  of  the  signals.  For  example,  the  interaction  of 
short  gravity  and  capillary  waves  with  surface  current  gradients,  such 
as  those  caused  by  Internal  waves  or  tidal  currents  flowing  over  bot¬ 
tom  topographic  features,  results  in  patterns  of  reflectivity  varia¬ 
tions  which  are  wavelength  dependent.  By  comparing  these  patterns  at 
different  wavelengths,  a  better  understanding  of  the  hydrodynamic 
Interactions  involved  can  be  gained,  and  these  wave-length  differences 
may  eventually  be  used  to  extract  more  detailed  oceanographic  informa¬ 
tion  from  the  images.  Full  exploitation  of  these  possibilities  re¬ 
quires  radiometric  calibration  of  the  SAR  data,  techniques  for  which 
are  currently  being  developed  [38]. 

A  second  type  of  mul ti -frequency  SAR  utilizes  the  interference 
between  the  surface  returns  at  two  closely  spaced  frequencies  in  order 
to  gain  information  about  the  surface  structure  on  a  spatial  scale 
which  is  related  to  the  difference  frequency  (Af) .  The  output  signal 
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from  such  a  SAR  system  is  proportional  to  the  surface  reflectivity 
spectral  density  at  the  surface  wavelength 


Xs  =  2(Af)  sin  9 


where  9  is  the  incidence  angle  and  c  is  the  speed  of  light.  The  prin¬ 
cipal  oceanographic  application  of  this  type  of  system  would  be  the 
measurement  of  spectral  perturbations  for  surface  wavelengths  not 
resolvable  by  conventional  SAR  systems.  Such  a  measurement  capability 
might  be  useful  for  the  detection  of  features  such  as  internal  waves 
which  cause  large  perturbations  of  waves  within  a  relatively  narrow 
spectral  region  due  to  resonant  interactions,  wave  trapping  or 
blocking. 

Preliminary  studies  have  indicated  the  feasibility  of  such  systems 
by  segmenting  the  range  bandwidth  of  a  conventional  SAR  system  into 
two  narrow  intervals  with  a  variable  spacing.  When  this  spacing  was 
set  to  correspond  to  the  wavelength  of  a  known  surface  feature,  ac¬ 
cording  to  Eq.  (9),  the  output  of  this  simulated  Af  system  was  ob¬ 
served  to  be  a  maximum.  The  range  of  observable  wavelengths  in  this 
demonstration  was  limited  by  the  signal  bandwidth  to  values  larger 
than  the  SAR  resolution.  However,  a  specially  designed  system  incor¬ 
porating  a  larger  frequency  difference  could  remove  this  restriction 
and  produce  more  useful  results. 

3.3  MULTIPLE  ANTENNA  SYSTEMS 

The  final  type  of  advanced  SAR  system  considered  in  this  report  is 
that  which  incorporates  multiple  antennas,  mounted  on  the  same  plat¬ 
form  but  displaced  in  either  the  along-track  or  across-track  direc¬ 
tions,  or  mounted  on  different  platforms.  The  purpose  of  these  ar¬ 
rangements  is  to  measure  special  attributes  of  the  scene  in  addition 
to  its  reflectivity,  such  as  the  surface  elevation  or  roughness,  or 
the  motion  of  the  surface  or  of  objects  within  the  scene. 
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A  two-antenna  SAR  system  for  measuring  target  motions  was  proposed 
by  Raney  [39]  using  a  phase  detection  method.  In  this  scheme,  coher¬ 
ent  pulses  are  transmitted  through  two  antennas  separated  in  the 
along-track  direction  by  a  distance  D,  with  the  pulse  fed  to  the 
second  antenna  delayed  by  an  amount 

St  =  D/V  (10) 

where  V  is  the  platform  velocity.  Thus,  the  two  pulses  are  transmit¬ 
ted  and  received  at  exactly  the  same  locations  relative  to  the  scene. 
If  the  scene  were  stationary,  therefore,  the  difference  between  the 
received  signals  would  be  zero,  whereas  if  a  given  scatterer  has  a 
radial  velocity  Vr,  the  received  signals  would  differ  by  a  phase 

bf  =  2kVr<5t  (11) 

where  k  is  the  radar  wavenumber,  and  the  magnitude  of  the  difference 
signal  would  be 

ffj  =  2  sin  (kVp«5t)  (12) 

where  I  si  is  the  magnitude  of  the  total  signal  received  from  the  scat¬ 
terer.  By  appropriately  choosing  the  antenna  separation  D,  any  given 
range  of  velocities  could  thus  be  measured.  For  example,  in  order  to 
measure  a  radial  velocity  on  the  order  of  0.1  m/s  from  an  earth¬ 
orbiting  platform  (V  ~  7  km/s)  with  an  X-band  system,  an  antenna  sepa¬ 
ration  on  the  order  of  35  meters  would  be  required. 

A  second  method  of  observing  the  motions  of  persistent  objects, 
such  as  sea  ice,  would  be  to  operate  a  squint-mode  SAR  with  two  anten¬ 
nas  having  an  angular  separation  bd,  so  that  the  same  area  is  imaged 
twice  with  a  time  difference 


At  =  (R/V)A0 


(13) 


where  R  Is  the  range  distance  and  V  is  the  platform  velocity.  If  the 
target  velocity  is  sufficient  to  cause  a  motion  through  several  reso¬ 
lutions  cells  during  this  time  interval,  the  velocity  would  be  measur¬ 
able  by  this  method.  For  a  squint-mode  SAR  operating  at  Seasat  alti¬ 
tudes  (R/V  ~  120  sec)  with  a  resolution  of  10  meters,  an  angular  sepa¬ 
ration  of  450  would  allow  velocities  on  the  order  of  0.5  m/s  to  be 
measured.  The  advantage  of  this  method  is  that  it  allows  both  compon¬ 
ents  of  the  velocity  to  be  measured  simultaneously. 

The  use  of  two  or  more  antennas  displaced  in  the  across-track 
direction  allows  the  possibility  of  measuring  surface  elevations  or 
roughness.  Basically,  this  is  because  two  or  more  instantaneous  mea¬ 
surements  of  range,  plus  a  measurement  of  the  along-track  coordinate, 
allows  each  scatterer  to  be  located  in  three  dimensions.  Several 
Implementations  of  this  concept  have  been  attempted,  Including  the 
hologram  radar  system  developed  at  ERIM  [40]  and  the  synthetic  inter¬ 
ferometer  radar  developed  at  Goodyear  Corp.  [41],  The  former  system 
used  an  array  of  100  receiving  antennas  arranged  in  the  across-  track 
direction,  the  signals  from  which  were  coherently  recorded  and  pro¬ 
cessed  to  yield  an  output  image.  Resolution  in  the  along-track  direc¬ 
tion  is  obtained  by  the  synthetic  aperture  technique.  By  transmitting 
two  closely-spaced  frequencies,  a  range  contouring  capability  was 
demonstrated  over  terrestrial  scenes.  A  modification  of  this  design 
could  conceivably  yield  useful  oceanographic  Information  such  as  wave 
heights  or  surface  roughness. 

The  Goodyear  system  [41]  uses  two  antennas  mounted  vertically  to 
obtain  two  independent  range  measurements,  and  also  utilizes  the  syn¬ 
thetic  aperture  technique  to  obtain  azimuth  resolution.  Summing  the 
signals  coherently  from  the  two  antennas  results  in  a  set  of  fringes 
or  nulls  at  fixed  depression  angles.  The  location  of  these  nulls  in 
the  output  Imagery  gives  information  on  both  the  range  and  the  depres¬ 
sion  angle,  which  allows  the  surface  elevation  to  be  computed  with  a 


known  ambiguity  factor.  This  concept  was  also  demonstrated  over  ter¬ 
restrial  scenes  and  may  conceivably  be  modified  or  generalized  to 
yield  oceanographic  information. 

The  term  bistatic  SAR  might  be  applied  to  any  of  the  multiple 
antenna  systems  described  above,  but  is  usually  reserved  for  the  case 
where  the  transmit  and  receive  antennas  are  located  on  different  plat¬ 
forms  and  the  antenna  separation  is  therefore  much  larger.  Previous 
discussions  of  bistatic  SAR  have  emphasized  the  possible  operational 
advantages  of  such  systems  for  real-time  and  continuous  monitoring 
applications  [42].  It  would  appear  that  additional  surface  elevation 
or  topographic  information  could  be  obtained  by  a  comparison  of  simul¬ 
taneous  monostatic  and  bistatic  images  of  the  same  area,  since  each 
contains  an  independent  measurement  of  the  range  distance  to  each 
target. 


4 

SUMMARY 


Strip-map  SAR  systems  are  useful  for  observing  a  variety  of 
oceanic  phenomena,  but  are  limited  primarily  by  their  sensitivity  to 
the  look  direction  of  the  SAR  relative  to  the  phenomenon  being  ob¬ 
served.  This  restriction  is  alleviated  by  the  use  of  angle  diversity 
SAR  techniques  which  allow  the  look  direction  to  be  varied  indepen¬ 
dently  of  the  direction  of  motion  of  the  SAR  platform.  These  tech¬ 
niques  may  also  prove  to  be  useful  for  gaining  information  on  scene 
motions  or  for  classifying  scene  elements  based  on  their  look  direc¬ 
tional  dependence.  Additionally,  angle  diversity  data  of  the  oceans 
can  be  processed  using  non-coherent  averaging  or  multiple-look  tech¬ 
niques  to  reduce  the  effects  of  background  noise.  This  background 
noise  is  multiplicative  in  nature  and  is  a  result  of  the  coherent 
nature  of  the  SAR.  Under  low  to  moderate  winds,  this  system  noise  or 
speckle  can  dominate  an  image  and  in  some  cases  mask  subtle  oceanog¬ 
raphy  features  of  interest.  _  This  is  a  particularly  acute  problem  in 
respect  to  detection  of  bottom  topography. 

Frequency  diversity  SAR  systems'  utilizing  multiple-wavelength 
transmissions  may  also  increase  the  accuracy  of  detection  and  classi¬ 
fication  of  surface  patterns,  based  on  the  wavelength  dependence  of 
the  surface  reflectivity.  Finally,  the  use  of  multiple  antenna  sys¬ 
tems  may  offer  the  possibility  of  measuring  surface  motions  at  high 
resolution,  or  of  providing  continuous  imaging  with  two  or  more  look 
di rections. 
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Numerical  Simulation  of  Synthetic  Aperture  Radar 
Image  Spectra  for  Ocean  Waves 

DAVID  R.  LYZENGA 


Abstract— A  numerical  model  for  predicting  the  synthetic  aperture 
radar  (SaR)  image  of  a  moving  ocean  surface  is  described,  and  results 
are  presented  for  two  SIR-B  data  sets  collected  otf  the  coast  of  Chile. 
Wave  height  spectra  measured  by  the  NASA  radar  ocean  wave  spec¬ 
trometer  (ROWS)  were  used  as  inputs  to  this  model,  and  results  are 
compared  with  actual  SiR-B  image  spectra  from  orbits  91  and  106. 
Additional  parametric  variations  are  presented  to  illustrate  the  effects 
of  nonlinearities  in  the  imaging  process. 

I.  Introduction 

HE  PROCESSES  involved  in  the  imaging  of  surface 
gravity  waves  by  synthetic  aperture  radar  have  been 
widely  discussed  and  debated  within  the  radio-oceano- 
graphic  community.  A  number  of  studies  have  been  ear¬ 
ned  out  to  investigate  these  processes  experimentally  [1]- 
[7]  as  well  as  theoretically  [8]-[  14|.  However,  because  of 
the  large  number  of  independent  variables  which  influ¬ 
ence  the  imaging  process,  the  experimental  difficulties  in¬ 
volved  in  making  the  appropriate  surface  measurements, 
and  differences  of  opinion  regarding  the  formulation  of 
the  theory,  a  clear  picture  of  the  limitations  of  SAR  for 
measunng  ocean  waves  has  been  slow  to  emerge. 

In  the  present  paper,  a  conceptually  simple  model  for 
the  SAR  imaging  process  is  discussed  and  the  numerical 
implementation  of  this  model  is  described.  Predictions  of 
the  SAR  image  spectra  for  various  environmental  condi¬ 
tions  and  SAR  system  parameters  are  presented,  and  com¬ 
parisons  are  made  with  measured  image  spectra  obtained 
during  the  SIR-B  experiment.  Additional  parametric  vari¬ 
ations  are  presented  to  illustrate  behavior  over  a  wider 
range  of  environmental  and  SAR  system  parameters,  and 
the  possibility  of  obtaining  wave  height  or  slope  spectra 
from  SAR  images  is  discussed. 

II.  Model  Description 

The  microwave  scattering  properties  of  the  ocean  sur¬ 
face.  including  the  phase  changes  due  to  surface  motions, 
may  be  described  in  terms  of  a  complex  reflectivity  p< x. 
. .  n  where  x  is  the  along-track  coordinate,  y  is  the  across- 
track  (or  range)  coordinate,  and  r  is  time.  The  signal  re¬ 
ceived  by  a  synthetic  aperture  radar  within  a  given  range 
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ceil  (assuming  the  signal  has  been  resolved  in  range,  and 
neglecting  range  resolution  effects)  may  then  be  written 
as 

st  y.  r)  =  \  pix.  y,  r)  a,(x  -  l Vt)  e~-'2knx-'-n  dx  (1) 

where  a ,(.r  -  Vt)  represents  the  antenna  gain  pattern,  k  is 
the  electromagnetic  wavenumber,  and  r(x.  y.  D  is  the  in¬ 
stantaneous  range  distance  from  the  scattering  element  at 
(.r.  y)  to  the  SAR  platform,  which  is  located  at  x  =  Vt.  y 
=  0,  and  c  =  A.  Assuming  that  the  antenna  pattern  limits 
the  illumination  to  an  along-track  distance  interval  that  is 
much  smaller  than  the  range  distance,  this  range  distance 
may  be  approximated  by 

_  (x  -  Vrr 

rix.  v,  t)  =  R  +  - — -  (2) 

where  R 1  =  y:  +  A2  The  integrand  in  (1).  which  may 
be  interpreted  as  the  contribution  of  a  surface  element  at 
(.r.  y)  to  the  total  received  signal,  can  then  be  written  as 

s(x.  y.  t)  =  p{x.  y.  t)  u,(.r  -  Vt)  e  """  ‘l '  '  R  < 3) 

neglecting  the  constant-phase  term  e'ZliR  that  is  common 
to  all  the  surface  elements  considered  in  the  integration. 

The  complex  image  amplitude  is  obtained  by  convolv¬ 
ing  the  recorded  signal  with  the  reference  function  lux'  I. 
i.e. 

i(x'.y')  =  j  s(y'.t)h(x'  -  Vt)  Jr  i4) 

where  the  reference  function  has  the  form 

lux'  -  Vt)  =  a,<.r'  -  Vt)  e  h"  ~‘r';  (5) 

in  which  uJx'  -  Vt)  is  the  amplitude  weighting  and  h  is 
the  quadratic  phase  factor  or  focus  parameter,  which  is 
normally  chosen  as 

b  =  k/R  ib) 

in  order  to  cancel  the  quadratic  time  dependence  in  the 
phase  of  the  recorded  signal.  If  this  is  so  chosen,  the  con¬ 
tribution  of  the  surface  element  at  (x.  y)  to  the  complex 
image  amplitude  at  (r  \  y’)  may  be  written  as 

Kx.x'.v')  =  p(.r.  y.  oj)  e  M*  R 
where  v  =  v’  i neglecting  range  resolution  effects)  and 

pt  v.  v.  oj)  =  |  pl.t.  v.  t)  <nx.  v  .  t)  c  ‘  Jr  iS) 


0106-2X02  86-  1  100-0X63S0I  00  |0X6  IEEE 


A- 3 


0  r 

v...- . . 


864 


IF.F.E  TRANSACTIONS  ON  GEOSCIENCE  AND  REMOTE  SENSING.  VOI  GF.  24.  NO  h.  NOVEMBER  IVX6 


with 


Scene 


SAR  Image 


a(.r,  x'.  t)  =  a\(x  -  Vt)  a2(x'  -  Vt)  (9) 

and 

w  =  2kV(x  -  x')/R.  (10) 

As  discussed  in  Hasselmann  et  al.  [14],  p(x.  v,  w)  may 
be  interpreted  as  the  finite-resolution  Fourier  transform  of 
the  surface  reflectivity,  with  the  resolution  being  deter¬ 
mined  by  the  width  of  the  weighting  function  a(x ,  r',  f). 
For  stationary  scenes,  p(. x,  y.  w)  is  centered  at  <*>  =  0  and 
has  a  width  inversely  proportional  to  the  width  of  a{x,  ,c'. 
t).  For  randomly  moving  surfaces,  p(x.  y.  w)  may  be  con¬ 
sidered  as  a  random  vanable  that  has  non-zero  values  over 
a  range  of  frequencies  determined  by  the  surface  motion 
and  the  radar  wavelength,  and  that  has  frequency  corre¬ 
lation  statistics  determined  by  a(x,  x  \  /). 

If  the  surface  reflectivity  is  assumed  to  be  spatially  un¬ 
correlated,  i.e. 


<p( x.  y.  t)  p*(: x\  y\  r')> 

=  <p(x,  y,  t)  p*(x'.  y’.  r')>  <5(.t  -  ,r')  <5(  v  -  v') 

(11) 

where  the  brackets  denote  ensemble  averaging  [14],  the 
expected  value  of  the  image  intensity  is  given  by 

</U'.  y')>  =  (  \  l(x-  x ’•  v')  dx  ^ 

=  j  fn(x.  y,  2 kV(x  -  x’)/R)  dx  (12) 


where  y  =  y '  and 
fD(x,  v,  oj)  = 


•> 

j  p(x.  y,  r)  u(x.  x',  t )  e  dt  ^ 


(13) 


is  the  finite-resolution  Doppler  spectrum  of  the  surface 
reflectivity. 

The  simulation  procedure  suggested  by  this  model  is 
therefore  to  compute  the  Doppler  spectrum  of  the  surface 
reflectivity  at  each  grid  cell  in  the  scene,  and  then  to  map 
the  radar  cross  section  onto  the  image  using  the  Doppler 
spectrum  to  control  the  mapping,  as  illustrated  in  Fig.  I . 
The  Doppler  spectrum  is  approximated  by  a  Gaussian-Iike 
function  (actually  a  polynomial  form  is  used)  with  a  cen¬ 
ter  location  and  a  width  The  total  area  under  the 
spectrum  is  equal  to  the  radar  cross  section  <r0U.  v). 

The  center  frequency  w0  and  width  a ^  may  be  expressed 
in  terms  of  the  surface  velocity  as 

w0  =  2  kVr  (14) 

and 

=  2  ka„  (15) 

where  V,  is  the  mean  radial  velocity  of  the  surface  and  a 
is  the  standard  deviation  of  the  radial  velocities  within  a 


Fig.  1.  Illustration  of  azimuthal  shifting  and  smearing  effects  due  to  sur¬ 
face  motions:  Xx  -  iR/V)  V,  and  hx  -  2 (R/V)o,..  where  V,  is  the  mean 
radial  velocity  and  <7,.  is  the  standard  deviation  of  radial  velocities  within 
a  given  resolution  cell. 


given  resolution  cell.  These  velocities,  and  the  radar  cross 
section,  are  computed  from  the  wave  height  spectrum  as 
follows. 

Values  of  the  wave  height  spectrum,  or  surface  eleva¬ 
tion  spectral  density,  are  assumed  to  be  specified  on  a  reg¬ 
ular  grid  of  spatial  wavenumbers 

k„,  =  —  (m  -  M/2).  m  =  1,  •  •  •  ,  M  -  1  (16) 

Lx 


and 


kn  =  —(n  -  N/2).  n  =  1,  •  •  •  ,  ,V  -  1  (17) 

L\ 

where  Lx  and  Lv  are  the  dimensions  of  the  scene  in  the 
along-track  and  across-track  directions,  respectively.  Each 
of  these  wavenumbers  is  assumed  to  correspond  to  a  freely 
propagating  deep-water  wave  of  amplitude  amn  and  fre¬ 
quency 

=  (gk)uz  (18) 

where  g  is  the  gravitational  acceleration  and  k  is  the  mag¬ 
nitude  of  the  wavenumber.  Thus,  the  surface  elevation  at 
time  r,  for  a  given  realization  of  the  spectrum,  may  be 
written  as 

If  —  i  v-  i 

Z{X.  y.  r)  =  S  S  a„„ 

m  =  I  n  =  I 

•  cos  (k,„x  +  k„v  -  u,nnr  *  o.„„)  ( 19) 

where  the  phases  <t>mn  are  randomly  selected  from  a  uni¬ 
form  distribution  over  (0.  2ir).  The  amplitudes  are 
proportional  to  the  square  root  of  the  spectral  density  at 
( km ,  kn),  with  the  normalization  chosen  to  produce  the  de¬ 
sired  value  for  the  significant  wave  height 

=  4<z:>'':  (20) 

where  (c)1  :  is  the  rms  surface  elevation. 

The  instantaneous  radial  velocity  of  the  surface  is  com¬ 
puted  as  the  sum  of  the  orbital  velocities  of  each  compo- 
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nent.  i.e. 

If-  i  v-  i 

Vr(Xt  V.  t)  *— '  Z  ^mn^mnSmn 

i  n  =  l 

•  cos  [kM.r  +  fc,  V  -  aj,„/  +  oTO  +  (21) 


The  radar  cross  section  is  calculated  in  terms  of  a  linear 
modulation  transfer  function,  as  discussed  by  Alpers  et 
al.  [  15].  Thus,  the  radar  cross  section  (averaged  over  A.t, 
Ay.  and  T)  is  given  by 

o0(x.  y,  t)  =  50(1  +  /(*.  y.  0|  (28) 


where 


gm  =  [ (k^Jky  sin:  9  +  cos2  0]‘ 


0V  =  tan' 


k  cos  9 
k.  sin  9 


and  9  is  the  incidence  angle. 

Assuming  that  scatterers  are  uniformly  distributed  over 
the  water  surface,  the  average  radial  velocity  over  each 
grid  cell  of  dimensions  A.r.  Ay  and  over  the  integration 
time  Tis  given  by  (21)  with  am„  replaced  by  amnbmn.  where 


The  radial  velocity  variance  is  calculated  as  the  sum  of 
the  contributions  from  each  resolved  wave  component  plus 
the  variance  due  to  the  subresolution-scale  waves,  i.e. 

j  *  -  i  v  -  i 

<j;(.r.  v.  t)  -  cr2,  +•  -  Z  Z  ainuii„gi„ 

^  *n  =*  I  n  ■  I 

[  ^  ’HPI 

•  cos2  (km. c  +  k,  V  -  <a>m„t  +  <t>mn  +  <j>,.)! 


where 


sin  (k„  A.r)  sin  (knSy)  sin  (wmn  T) 
kmAx  k„  Av  w„„7’ 


o:t  =  1  1  u)'[sin2  o  sin2  9  +  cos2  9]  S(k.  <t>)k  dk  J<t> 

VT\it  O  f) 

(27) 

where  S(k.  ©)  is  the  spectrum  of  the  subresoiution-scale 
waves.  9  is  the  incidence  angle,  and  kr  is  the  electromag¬ 
netic  cutotf  wavenumber  in  the  two-scale  model  [14],  The 
summation  in  (25)  incorporates  both  spatial  and  temporal 
variations  in  the  radial  velocity  due  to  long  waves.  In  the 
limit  as  k  —  0  and  T  —  0,  the  contribution  of  a  given 
wave  component  to  the  velocity  variance  is  equivalent  to 
that  calculated  for  a  surface  with  a  constant  acceleration. 
However,  these  contributions  are  added  incoherently,  so 
the  spatial  variations  in  a2,  calculated  from  (25)  are  not 
necessarily  the  same  as  those  obtained  from  a  calculation 
of  the  acceleration  effect. 


where  <j0  is  the  mean  radar  cross  section  for  the  entire 
scene,  and 

v- i  v- i 

fix.  v.  r)  =  Z  Z  a^b^R^ 

m  ~  I  n  *  | 

•  cos  [k„.r  +  kny  -  +  <&„„  -  oR\ 

(29) 

where  R„„  and  are  the  modulus  and  phase,  respec¬ 
tively,  of  the  radar  modulation  transfer  function.  This 
transfer  function  includes  the  effects  of  changes  in  the  lo¬ 
cal  slope  of  the  surface  (tilt  modulation)  as  well  as  changes 
in  the  small-scale  surface  roughness  (hydrodynamic  mod¬ 
ulation).  Note  that  this  modulation  transfer  function  is 
equivalent  to  that  discussed  by  Alpers  et  al.  [15]  and  is  k 
times  the  dimensionless  MTF  defined  by  Plant  et  al.  [16]. 
For  the  case  of  tilt  modulation,  this  transfer  function  may 
be  written  as 

=  k„oo(W°o(0)  (30) 

and 

9>r  —  ±  ~  (31) 


where  ao(0)  is  the  radar  cross  section  of  the  surface  as  a 
function  of  incidence  angle,  and  the  sign  of  <t>R  is  positive 
for  look  directions  in  the  positive  y-direction  and  negative 
for  look  directions  in  the  negative  y-direction  [15]  Using 
a  Bragg  scattering  model  for  ao(0 )  and  assuming  a  per¬ 
fectly  conducting  surface  and  a  Phillip's  spectrum  for  the 
short  waves,  the  tilt  modulation  transfer  function  may  be 
written  as 


4  cot  9 
l  ±  sin2  9 


(32) 


where  the  (  +■ )  plus  sign  holds  for  vertical  polarization  and 
the  (-)  minus  sign  is  for  horizontal  polarization.  Note 
that  this  transfer  function  goes  to  zero  for  azimuth-trav¬ 
eling  waves  (k,  =  0). 

The  hydrodynamic  contribution  to  the  modulation 
transfer  function  involves  several  possible  mechanisms, 
including  straining  of  short  waves  by  long-wave  orbital 
motions,  modulation  of  the  air  flow  across  the  long  waves, 
and  nonlinear  effects  such  as  wave  breaking  Not  all  of 
these  mechanisms  are  well  understood  on  a  theoretical  ba¬ 
sis.  but  intuitively  one  might  expect  the  effects  to  depend 
on  environmental  conditions.  Measurements  of  the  mod¬ 
ulation  transfer  function  indicate  a  wind-speed  depen 
dence  at  .Y-band  radar  frequencies,  but  no  apparent  wind- 
speed  dependence  at  T-band  [16] 
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The  simulation  procedure  described  above  yields  the 
ensemble-averaged  or  expected  value  of  the  image  inten¬ 
sity,  which  will  be  denoted  by  /0U.  v)  in  the  following. 
The  image  intensity  /(. t,  v)  observed  with  an  actual  SAR 
is  influenced  by  coherent  interference  or  speckle  effects, 
which  may  be  modeled  as  a  multiplicative  noise  process, 
i  e. 


l(x,  v)  =  /0(.r.  v)[l  +  e(.x.  y)J 


(33) 


where  et.t.  y)  is  a  zero-mean  random  variable  which  has 
a  variance  of  1  ,V;,  where  .V,  is  the  number  of  “looks” 
used  in  processing  the  data  [17],  e(.t.  v)  also  has  spatial 
correlation  statistics  that  are  dependent  on  the  SAR  sys¬ 
tem  bandwidth  and  are  independent  of  scene  motions  [18]. 

If  H.x.  >•)  is  computed  on  a  spatial  grid  with  spacings 
larger  than  or  equal  to  the  nominal  SAR  resolution,  the 
values  of  et.t.  y)  may  be  considered  to  be  uncorrelated 
from  point  to  point.  Thus,  the  spectrum  of  <s(.c.  v).  which 
will  be  denoted  by  S,{kx,  ks).  is  white.  The  spectrum  of 
H.x.  v)  can  then  be  written  as 


S,{k,.  kv)  =  50(£r,  k,)  +  Sn(k„  fcv) 


(34) 


where  S^kx.  kv)  is  the  spectrum  of  /0(.t.  y),  and  S„(kx.  kx) 
is  the  convolution  of  S0(kx.  k v)  with  S,(kx,  kx).  which  re¬ 
sults  in  another  white  spectrum  [19].  Thus,  the  expected 
value  of  the  image  spectrum  is  equal  to  the  spectrum  of 
/,)(. t.  y>  plus  a  constant  noise  floor.  The  mean  value  of  this 
noise  floor  can  be  calculated  using  Parseval’s  theorem  and 
the  fact  that 


l/.V,) 


(35) 


</->  =  </r,X 

which  follows  readily  from  (33).  where  the  triangular 
brackets  refer  to  average  values  over  the  entire  image. 
From  Parseval's  theorem,  the  average  overall  wavenum¬ 
bers  of  each  term  in  (34)  is  proportional  to  the  corre¬ 
sponding  term  in  (35)  Thus,  the  average  or  expected 
value  of  the  noise  floor  is 


<S„>  -  <S„>/M 


(36) 


where  <  50>  is  the  average  value  of  Sq(kx,  kj  over  all 
wavenumbers,  including  k,  =  0.  kr  =  0. 

A  measure  of  the  wave  contrast  or  detectability  is  given 
by  the  peak-to-background  ratio,  defined  as  the  peak 
spectral  density  in  the  image  spectrum  divided  by  the  noise 
spectral  density.  From  the  derivation  presented  above,  this 
ratio  may  be  expressed  as 


PBR  = 


<S„> 


(S<))n 


<S0> 


.V, 


(37) 


where  is  the  peak  value  in  the  spectrum  of  v). 

A  comparison  of  predicted  values  of  this  parameter  with 
observed  values  is  presented,  along  with  comparisons  of 
predicted  and  observed  image  spectra,  in  the  following 
section. 

III.  Model  Results 

The  simulation  model  described  above  has  been  exer¬ 
cised  using  wave  height  spectra  measured  by  the  NASA 


Radar  Ocean  Wave  Spectrometer  [20]  during  the  SIR-B 
Southern  Oceans  experiment  [7].  The  results  of  these  cal¬ 
culations  are  presented  here  and  are  compared  with  actual 
SIR-B  image  spectra.  Additional  results  are  also  pre¬ 
sented  illustrating  the  imaging  behavior  for  parameter 
variations  not  encountered  during  the  SIR-B  experiment. 

Data  were  collected  by  the  SIR-B  sensor  during  orbit 
91  on  October  1 1 .  1984  off  the  coast  of  Chile,  and  wave 
height  spectra  were  measured  nearly  simultaneously  in  the 
same  area  by  the  ROWS  system  aboard  a  NASA  P-3  air¬ 
craft.  Detailed  descriptions  and  intercomparisons  of  these 
data  sets  are  given  in  Beal  et  al.  [7],  Contour  plots  of  the 
ROWS  wave  height  and  wave  slope  spectra  are  shown  in 
Fig.  2(a)  and  (b).  Note  that  the  coordinate  system  is  re¬ 
ferred  to  the  SIR-B  ground  track  direction,  which  is  hor¬ 
izontal  on  these  plots. 

The  wave  height  spectrum  shown  in  Fig.  2(a)  was  used 
in  the  SAR  model  described  in  the  previous  section,  along 
with  the  other  parameters  listed  in  Table  I.  The  radar 
modulation  transfer  function  was  evaluated  from  (32)  at  d 
=  25°,  for  horizontal  polarization.  This  yields  a  dimen¬ 
sionless  MTF  having  a  magnitude  of  10  and  a  phase  of 
90°.  The  contribution  to  the  radial  velocity  variance  from 
subresolution-scale  waves  was  computed  using  a  Pierson- 
Moskowitz  spectrum  [21]  with  a  wind  speed  of  12  m/s. 
The  velocity  variance  resulting  from  this  calculation  is 
approximately  0.21  (m/s);.  The  velocity  variance  contrib¬ 
uted  by  waves  longer  than  50  m  is  a  function  of  position, 
but  has  an  average  value  of  0. 19  (m/s):  tor  this  case.  The 
spectrum  of  the  simulated  SAR  image  intensity  (/„)  is 
shown  in  Fig.  2(c),  and  the  actual  SIR-B  image  spectrum 
is  shown  in  Fig.  2(d). 

An  intercomparison  of  the  spectra  shown  in  Fig  2 
shows  that  the  slope  spectrum,  the  simulated  image  spec¬ 
trum,  and  the  actual  image  spectrum  are  all  similar  in 
shape.  The  similarity  between  the  simulated  image  spec¬ 
trum  and  the  slope  spectrum  results  from  the  use  of  a  di¬ 
mensionless  modulation  transfer  function  that  is  indepen¬ 
dent  of  wavelength.  Since  both  the  dominant  and 
secondary  wave  systems  are  nearly  range-traveling,  sur¬ 
face  motion  effects  are  minimal  and  the  image  is  domi¬ 
nated  by  real  modulation  effects.  The  fact  that  the  actual 
image  spectrum  has  a  similar  shape  supports  the  assump¬ 
tion  of  a  wavelength-independent  MTF  Although  there 
are  slight  differences  in  the  locations  and  relative  magni¬ 
tudes  of  the  primary  and  secondary  peaks,  these  differ¬ 
ences  are  felt  to  be  within  the  limits  of  experimental  ac¬ 
curacy  . 

The  ratio  of  the  peak  spectral  density  to  the  mean  spec¬ 
tral  density  for  the  simulated  image  is 


(■5o)ma\ 

<5„> 


=  9.5 


(38) 


for  this  case.  From  (37),  this  implies  a  peak-to-back¬ 
ground  ratio  of  10.5  (for  ,V,  =  1)  when  speckle  effects  are 
included.  For  the  actual  SIR-B  image  spectrum,  a  com¬ 
parable  value  was  obtained  by  computing  the  background 
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Fig.  2.  Comparison  of  (al  wave  height.  (b>  wave  slope.  (cl  simulated  im¬ 
age.  and  Id)  actual  SIR-B  image  spectra  for  orbit  91  lOctober  11. 
19841, 


TABLE  I 

Environmental  and  SAR  System  Parameters  for  Orbit  91 


Significant  <«  ave  -feignt 

5.0  ■n 

.aci!  rttrui  "oeea 

'2  "/S 

SI3-3  -‘e-iti  -Ig 

54.2 *r 

ir.'ienr.8  ingle 

25.1" 

3 iatf irm  A  I  *  1 C jde 

23 7 

3  ’  3  /<*  Iqc  •  \f 

’  =  00  n/s 

!nteg*“*t-on  ’une 

3. 56  sec 

3,<e>  Soac*ng 

25  <  25  m 

spectral  density  along  the  right-hand  side  of  the  spectrum 
(at  fc,  =  2 ir / 100  m).  The  ratio  of  the  peak  spectral  density 
to  this  background  value  is  8.8. 

The  second  SIR-B  data  set  for  which  comparisons  have 
been  made  was  collected  during  orbit  106  on  October  12. 
1984.  Contour  plots  for  the  ROWS  wave  height  and  slope 


spectra  for  this  pass  are  shown  in  Fig.  3(a)  and  (b).  The 
parameters  used  in  the  SAR  model  for  this  pass  are  given 
in  Table  II.  The  same  value  was  assumed  for  the  radar 
MTF  as  in  the  previous  case.  However,  this  parameter 
has  almost  no  effect  on  the  results  of  the  simulation,  since 
the  dominant  waves  are  traveling  in  the  azimuth  direction 
where  the  tilt  MTF  goes  to  zero.  The  radial  velocity  vari¬ 
ance  due  to  subresolution-scale  waves  was  assumed  to  be 
the  same  as  in  previous  case,  and  the  variance  contributed 
by  longer  waves  is  approximately  0.09  (m/s):  for  this 
case.  The  simulated  and  actual  SIR-B  image  spectra  are 
shown  in  Fig.  3(c)  and  (d),  respectively. 

The  simulated  and  actual  image  spectra  are  again  quite 
similar  in  appearance  to  the  slope  spectrum  for  this  data 
set.  However,  the  similarity  between  the  simulated  image 
spectrum  and  the  slope  spectrum  arises  from  a  completely 
different  mechanism  in  this  case.  The  modulation  of  the 
backscattered  power,  as  calculated  from  surface  tilt  ef- 
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Fig.  3.  Comparison  of  (a)  wave  height.  (b)  wave  slope,  (c)  simulated  im¬ 
age.  and  (d)  actual  SIR-B  image  spectra  for  orbit  106  (October  12. 
1984). 


TABLE  II 

Environmental  and  SAR  System  Parameters  for  Orbit  106 


Significant  Wave  Height 

3.5  i» 

Local  Aina  Soeeo 

1?  n/i 

SIR-3  Heading 

95. VT 

Incidence  Angle 

18.3' 

3 latfonn  Altitude 

237  <m 

3Utfor»n  Velocity 

7500  'n/s 

Integration  'ime 

0.63  sec 

3,*el  Spacing 

25  <  25  * 

fects.  is  nearly  zero  for  the  dominant  azimuth-traveling 
waves  present  in  this  pass.  Thus,  the  calculated  image 
intensity  modulation  is  due  almost  entirely  to  surface  mo¬ 
tion  (velocity  bunching)  effects.  The  apparent  agreement 
between  the  calculated  and  observed  spectra  provides 
strong  support  for  the  assumed  velocity  bunching  mech¬ 
anism,  although  it  does  not  rule  out  the  possibility  of  other 


modulation  mechanisms  such  as  nonlinear  hydrodynamic 
effects  which  might  produce  backscatter  modulations  for 
azimuth  traveling  waves. 

The  simulated  image  spectrum  for  this  case  has  a  nor¬ 
malized  peak  spectral  density  of 


which  yields  a  peak-to-background  ratio  of  4. 1  from  (37V 
The  peak-to-background  ratio  of  the  actual  image  spec¬ 
trum,  calculated  using  the  same  procedure  as  described 
above,  is  2.8.  The  difference  between  the  observed  and 
predicted  values  of  the  PBR  may  be  due  to  additional  noise 
sources  affecting  the  background  spectral  density.  The 
rather  large  difference  between  the  PBRs  for  the  two  cases 
considered  (8.8  for  orbit  91  versus  2.8  for  orbit  106)  is 
predicted  by  the  simulation  model.  This  qualitative  agree¬ 
ment,  along  with  the  similarity  in  shape  between  the  sim- 
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ulated  and  actual  image  spectra,  lends  confidence  to  the 
simulation  procedure.  More  extensive  comparisons  with 
observations  over  a  wider  range  of  conditions  are  needed 
to  fully  validate  the  model,  however. 

The  effects  of  surface  motions  are  dependent  on  the 
SAR  system  parameters  (primarily  the  R/V  ratio)  as  well 
as  the  surface  conditions  (as  described  by  the  wave  height 
spectrum)  and  are  most  critical  for  waves  traveling  in  the 
azimuth  direction.  In  order  to  illustrate  these  effects  over 
a  wider  range  of  conditions  than  occurred  during  the  SIR- 
B  experiment,  the  simulation  model  was  operated  for 
R/V  values  of  30.  60.  and  120  s  and  for  wave  spectra 
having  peak  wavelengths  of  100.  200.  and  400  m  in  the 
azimuth  direction.  The  pixel  spacing  was  chosen  to  be  25 
x  25  m  in  each  case,  and  the  integration  times  were  as¬ 
sumed  to  be  0.56,  1.13,  and  2.26  s  for  the  three  R/V  val¬ 
ues  selected.  These  parameters  are  typical  of  spacebome 
L-band  SAR  systems  such  as  SIR-B  (with  R/V  ranging 
from  30  to  60  s)  and  Seasat  or  ERS-1  (with  R/V  =  120 
s).  The  incidence  angle  was  taken  as  25°  in  each  case, 
and  the  radar  modulation  transfer  function  was  assumed 
to  be  zero. 

For  the  purpose  of  these  illustrative  calculations,  a  sim¬ 
ple  form  was  assumed  for  the  wave  height  spectrum,  as 
given  by 


A  exp  [  -2(k,y  k  cos  o):] 
[k:  -  2kk}  cos  ®  +  ton'll2 


ir 

-  <  <t>  < 


x 


(40) 


This  spectrum  has  a  peak  spectral  density  at  k  =  kn,  <b  = 
0.  and  has  a  frequency  spectrum  nearly  identical  to  that 
proposed  by  Pierson  and  Moskowitz  [21].  if  A  =  0.0013, 
k0  =  0.59  g/u~  and  kt  =  0.4  kn.  This  spectrum  also  has 
a  frequency-dependent  angular  spreading  which  is  similar 
to  that  inferred  by  Hasselmann  et  al.  [22]  from  the  JON- 
SWAP  1973  data  for  frequencies  above  the  spectral  peak, 
and  approaches  the  semi-isotropic  Phillips  spectrum  fori- 
»  /c0.  Thus,  this  form  appears  to  be  a  reasonable  rep¬ 
resentation  of  the  wave  spectrum  for  a  fully  developed 
sea. 

The  wave  spectrum  defined  by  (40).  with  parameters 
selected  to  match  the  Pierson-Moskowitz  spectrum,  im¬ 
plies  a  significant  wave  height  ( H,)  of  approximately  2 
percent  of  the  peak  wavelength  (which,  incidentally,  is 
roughly  30  percent  larger  than  2xg/u>;.  where  ujp  is  the 
peak  frequency).  In  order  to  investigate  the  effects  of 
changes  in  wave  slope,  runs  were  also  made  with  the  value 
of  A  reduced  to  yield  a  ratio  Hv  \p  =  0.01.  where  = 
2x/k0  is  the  peak  wavelength.  The  wave  spectrum  used 
for  the  subresolution-scale  waves  was  also  reduced  by 
one-half  for  these  runs. 

The  results  of  these  simulations  are  shown  in  Figs.  4 
through  6.  Fig.  4  shows  the  wave  height  and  slope  spectra 
and  the  corresponding  image  spectra  for  the  case  = 
200  m  and  R/V  =  30  s.  Fig.  5  shows  the  ratio  of  the  image 
spectrum  to  the  wave  height  spectrum  (hereafter  referred 
to  as  the  SAR  transfer  function)  along  the  azimuth  wavc- 


number  axis  for  the  same  R/V  ratio.  The  straight  line  in 
this  plot  indicates  the  velocity  bunching  transfer  function 

i  R\~ 

Ri  =  (  —  j  gk}  cos:0  G:<0.  o)  (41) 

defined  by  Alpers  et  at.  [15],  where 

Gl(6.  0)  =  sin:  9  sin'  0  *  cos'  9  (42) 

and  0  =  0  for  the  case  under  consideration.  The  SAR 
transfer  function  obtained  from  the  simulations  ap¬ 
proaches  this  line  (which  indicates  the  mapping  is  in  the 
“linear"  regime)  when  the  significant  wave  height  is  suf¬ 
ficiently  small.  For  these  cases  (with  R/V  =  30  sec  and 
\p  <  400  m).  the  modulation  produced  by  velocity 
bunching  is  larger  than  that  which  would  be  produced  by 
tilt  modulation  if  the  waves  were  range-traveling.  In  the 
cases  where  the  slope  of  the  transfer  function  is  near  2. 
the  SAR  image  spectrum  resembles  the  wave  slope  spec¬ 
trum.  as  shown  in  Fig.  4(d). 

When  the  wave  height  or  R/V  ratio  increases,  the  trans¬ 
fer  function  begins  to  deviate  from  the  value  indicated  by 
(41).  For  these  cases,  the  transfer  function  can  be  approx¬ 
imately  modeled  as  the  product  of  (41)  with  an  “azimuth 
falloff"  function  that  is  presumed  to  be  related  to  the  sur¬ 
face  velocity  distribution  [6],  (23).  Over  certain  wave- 
number  regions,  the  product  of  these  two  functions  may 
be  roughly  constant,  as  in  the  (100,  1)  and  (200.  4)  cases 
shown  in  Fig.  5.  For  this  subset  of  cases,  the  SAR  image 
spectrum  more  closely  resembles  the  wave  height  spec¬ 
trum.  as  shown  in  Fig.  4(c).  When  the  wave  height  or  R/ 
V  ratio  increases  further,  the  transfer  function  is  domi¬ 
nated  by  the  azimuth  falloff  effect,  and  the  image  spec¬ 
trum  becomes  highly  distorted. 

An  example  of  the  effects  of  changing  R/V  are  shown 
in  Fig.  6.  This  plot  of  the  SAR  transfer  function  for  = 
400  m  and  H,  =  4  m  indicates  a  gradual  transition  from 
the  linear  regime  (for  R/V  =  30  s)  to  a  moderately  non¬ 
linear  regime  (for  R/V  -  120  s).  The  azimuth  falloff  func¬ 
tion  for  these  cases  can  be  approximated  by 

p(k,)  =  exp  (  - a]k ;)  (43) 

where  a,  has  the  value  17-m  for  R/V  =  30  s.  34  m  for 
R/V  =  60  s.  and  64  m  for  R/V  =  120  s.  These  values  are 
consistent  with  a  velocity  variance  of  approximately  0.32 
(m/s):,  assuming  the  effective  azimuthal  impulse  response 
function  has  the  form 

/  I  Vx  \ 

>'■"  ’  “P  ( lb.)  ,M' 

where  a:,  is  the  velocity  variance,  as  discussed  in  Section 
II.  This  value  is  more  than  twice  as  large  as  the  velocity 
variance  calculated  from  (25)  for  this  case,  but  is  approx¬ 
imately  equal  to  the  total  velocity  variance  for  the  as¬ 
sumed  wave  height  spectrum.  Thus,  it  appears  that  sur¬ 
face  motions  associated  with  all  wavelength  scales 
influence  the  azimuth  falloff  effect.  This  is  consistent  with 
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Fig.  4.  Companson  of  (a)  wave  height,  (b)  wave  slope,  and  (c)  simulated 
image  spectra  for  V  =*  30  s,  X„  *  200  m.  and  W,  =  4  m.  and  (dl 
simulated  image  spectrum  for  R/V  *  30  s.  \r  =  200  m.  and  H,  =  2  m. 


rig  5  Ratio  of  image  spectrum  to  wave  height  spectrum  along  azimuth 
wavenumber  axis  for  R/V  *  30  s.  Labels  indicate  dominant  wave  length 
\XP)  and  significant  wave  height  (tft)  Solid  line  indicates  the  velocitv 
hunching  transfer  function 


Fig  5  Ratio  oi  image  spectrum  to  wave  height  spectrum  along  azimuth 
wavenumber  axis  tor  R  V  =  30.  bO.  and  120  s.  with  *  400  m  and  H 
=  4  m  Solid  curves  indicate  values  obtained  by  simulations.  Dotted 
curves  indicate  analytical  fits  discussed  in  text 
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observations  by  Beal  et  al.  [6]  indicating  that  the  azimuth 
falloff  is  dependent  on  the  significant  wave  height. 

An  azimuth  cutoff  wavelength  can  be  somewhat  arbi¬ 
trarily  defined  as 

K  =  *ai  (45) 

which  corresponds  to  the  wavelength  at  which  the  falloff 
function  reaches  a  value  of  approximately  0.018.  \c  is 
also  roughly  half  the  wavelength  at  which  the  total  SAR 
transfer  function  reaches  its  peak  value.  The  numerical 
values  of  \r  for  the  examples  shown  in  Fig.  6  are  approx¬ 
imately  equal  to  those  predicted  by  Beal  et  al.  [6]  for  a 
significant  wave  height  of  4  m.  Note  that  the  peak  value 
of  the  SAR  transfer  function  decreases  with  increasing 
R/V  for  a  given  wave  spectrum.  Thus,  for  wavelengths 
near  this  peak,  the  wave  contrast  actually  decreases  with 
an  increase  in  R/V.  while  the  reverse  is  true  for  much 
longer  wavelengths. 

IV.  Summary  and  Conclusions 

The  numerical  model  described  in  this  paper  indicates 
that  for  sufficiently  small  waveheight/waveiength  and 
R/V  ratios,  the  imaging  of  ocean  waves  by  synthetic  ap¬ 
erture  radar  may  be  considered  to  be  a  linear  process,  in 
the  sense  that  the  image  spectral  density  at  any  given 
wavenumber  (after  correcting  for  image  speckle  effects) 
is  proportional  to  the  wave  height  spectral  density  at  that 
wavenumber.  Assuming  that  the  constant  of  proportion¬ 
ality  (or  SAR  transfer  function)  can  be  estimated,  the  wave 
height  or  slope  spectrum  can  then  be  inferred  from  a  mea¬ 
surement  of  the  SAR  image  spectrum.  For  low-altitude 
spacecraft  sensors  such  as  SIR-B.  this  linear  regime  en¬ 
compasses  a  reasonably  large  range  of  wave  conditions. 
Predictions  of  image  spectra  for  two  SIR-B  data  sets  agree 
reasonably  well  with  actual  image  spectra  and  indicate  that 
the  imaging  is  within  the  linear  regime  for  these  cases. 

As  the  wave  height  or  R/V  ratio  increases,  the  imaging 
process  enters  a  regime  in  which  the  ratio  of  the  image 
spectrum  to  the  wave  height  spectrum  is  dependent  on 
some  integrated  property  of  the  wave  spectrum  itself  (pos¬ 
sibly  the  significant  wave  height).  In  this  case,  estimating 
the  wave  height  or  slope  spectrum  from  the  SAR  image  is 
obviously  more  difficult,  although  it  is  conceivable  that 
under  certain  circumstances,  an  iterative  procedure  might 
be  used  to  accomplish  this  inversion.  As  the  wave  height 
or  R/V  ratio  increases  further,  surface  motion  effects  cause 
increasing  distortions  in  the  image  spectrum  and  the  pos¬ 
sibility  of  estimating  the  wave  height  spectrum  from  the 
SAR  image  becomes  more  remote. 
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On  the  Estimation  of  Wave  Slope-  and 
Height- Variance  Spectra  from  SAR  Imagery 

FRANK  M.  MONALDO  and  DAVID  R.  LYZENGA 


Abstract—  A  procedure  is  described  for  using  synthetic  aperture  ra¬ 
dar  (SAR)  imagery  to  estimate  two-dimensional  ocean  wave  slope-  and 
height-variance  spectra.  The  logic  underpinning  the  procedure  is  based 
both  on  the  results  of  the  numerical  simulation  of  SAR  wave  imagery 
and  analytic  descriptions  of  the  SAR  imaging  process.  The  procedure, 
when  applied  to  SAR  imagery  of  waves  acquired  during  the  recent 
Shuttle  Imaging  Radar  Mission  (SIR-B),  is  shown  to  produce  spectra 
that  agree  with  independent  measures  of  both  the  two-dimensional 
slope-  and  height-variance  spectra.  The  implications  uf  these  results 
for  future  SAR  missions  aimed  at  measuring  ocean  waves  are  consid¬ 
ered. 

I.  Introduction 

HIGH-RESOLUTION  SAR  imagery  has  shown  the 
potential  of  measuring  ocean  wave  spectra.  Using 
data  from  both  spacebome  [I],  [2]  and  a  number  of  air¬ 
craft  [3],  (4)  platforms,  various  investigators  have  dem¬ 
onstrated  that  spectral  parameters  such  as  wave  number 
and  propagation  direction  extracted  from  SAR  image 
spectra  are  generally  consistent  with  independent  mea¬ 
sures  of  the  ocean  surface.  In  high  sea  states  though,  a 
rotation  of  the  wave  vector  toward  the  range  (cross-track) 
direction  can  occur  [1],  [2],  [5],  [6], 

In  spite  of  voluminous  imagery  from  both  the  Seasat 
SAR  and  various  aircraft  SAR  systems,  a  generally  ac¬ 
ceptable  procedure  for  generating  estimates  of  the  ocean 
wave  slope-  or  height-variance  spectrum  from  SAR  im¬ 
agery  has  not  been  found. 

The  basic  problem  is  in  specifying  the  relationship  be¬ 
tween  SAR  image  intensity  and  either  ocean  surface  wave 
height  or  slope.  Although  tower-based  measurements  of 
the  relationship  between  local  surface  slope  and  radar 
cross  section  have  yielded  valuable  information,  the  map¬ 
ping  of  surface  radar  cross  section  into  SAR  image  inten¬ 
sity  is  not  direct.  Because  a  SAR  relies  on  Doppler  infor¬ 
mation  to  achieve  high  azimuth  (along-track)  resolution, 
ocean  surface  movement  causes  the  radar  cross  section 
associated  with  a  particular  spot  on  the  ocean  surface  to 
be  mapped  into  image  intensity  at  a  variety  of  locations 
in  the  final  SAR  image. 

A  second  difficulty  associated  with  prescribing  a  pro¬ 
cedure  for  extracting  wave  spectra  from  SAR  imagery 
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is  the  lack  of  independently  measured  two-dimensional 
ocean  wave  spectra  over  scales  comparable  to  SAR  image 
spectra.  Without  such  independent  spectra,  it  is  impossi¬ 
ble  to  verify  ocean  slope-  or  height-variance  spectra  es¬ 
timates  from  SAR  imagery. 

These  two  difficulties  have  been  alleviated  by  two  re¬ 
cent  developments:  1)  a  sufficiently  complete  SAR  im¬ 
aging  theory  to  allow  complete  two-dimensional  model¬ 
ing,  both  analytic  [7]  and  numerical  [8],  [9],  of  the  SAR 
wave  imaging  process,  and  2)  new  methods  of  measuring 
the  two-dimensional  ocean  wave  spectrum  with  high 
spectral  resolution. 

Alpers  et  al.  [7]  have  described  a  linear  systems  ap¬ 
proach  to  the  SAR  imaging  process  that  provides  a  basis 
for  a  procedure  to  estimate  wave  height  spectra  from  SAR 
image  spectra.  Monaido  [6]  has  demonstrated  the  addi¬ 
tional  value  of  treating  the  indirect  and  formally  nonlinear 
mapping  of  surface  radar  cross  section  into  various  pixel 
locations  in  the  SAR  image  as  a  low-pass  filter  in  a  linear 
system  so  as  to  retain  the  advantages  of  tractability  inher¬ 
ent  in  a  linear  systems  approach. 

The  limits  to  the  effectiveness  of  a  linear  systems  ap¬ 
proach  to  SAR  wave  imaging  can  be  explored  through  the 
use  of  numerical  simulations.  SAR  wave  image  simula¬ 
tions  have  been  implemented  by  Alpers  [8]  and  Lyzenga 
[9].  The  simulations  begin  with  an  arbitrary  input  ocean 
spectrum,  construct  the  corresponding  ocean  surface 
height  topography,  use  models  for  the  modulation  of  the 
backscatter  as  function  of  wave  height,  and  include  the 
indirect  mapping  of  intensity  caused  by  surface  motion  to 
form  a  SAR  image.  A  spectrum  of  the  image  is  computed 
and  then  compared  to  the  input  ocean  spectrum.  In  this 
paper,  the  numerical  SAR  wave  image  modeling,  as  im¬ 
plemented  by  Lyzenga,  is  used  to  address  the  limits  of 
treating  the  effect  of  ocean  surface  motion  as  a  low-pass 
filter. 

In  the  seven  years  since  the  launch  and  subsequent  de¬ 
mise  of  Seasat.  at  least  two  airborne  techniques  for  the 
measurement  of  directional  ocean  wave  spectra  have  ma¬ 
tured:  the  surface  contour  radar  (SCR)  [10]  and  the  radar 
ocean  wave  spectrometer  (ROWS)  [11].  These  instru¬ 
ments  have  been  shown  to  reliably  measure  the  two-di¬ 
mensional  ocean  surface  wave  spectrum  with  wavenum¬ 
ber  and  angular  resolution  comparable  to  that  obtained 
from  a  typical  SAR  image  spectrum. 

The  Shuttle  Imaging  Radar  Mission  (SIR-B)  aboard  the 
shuttle  Challenger  permitted  the  execution  of  an  experi- 
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ment  where  SAR  imagery  was  obtained  in  regions  spa¬ 
tially  and  temporally  coincident  with  measurements  of  the 
ocean  wave  spectrum  by  the  SCR  and  ROWS  mounted 
aboard  a  NASA  P-3  aircraft.  Preliminary  results  from  the 
experiment  are  reported  in  [12]. 

The  P-3  also  was  equipped  with  two  other  wave  mea¬ 
suring  instruments:  1)  the  advanced  airborne  flight  exper¬ 
iment  (AAFE)  altimeter,  which  used  the  spreading  of  a 
narrow  nadir-transmitted  radar  pulse  to  estimate  the  total 
height  variance  of  the  surface,  and  an  airborne  optical  li- 
dar  (AOL),  which  provides  wave  height  profiles  along  the 
aircraft  flight  direction  that  can  be  subsequently  processed 
to  yield  one-dimensional  ocean  wave  spectra. 

In  October,  the  highest  sea  states  are  generally  found 
at  far  southern  latitudes.  For  this  reason,  the  experiment 
was  staged  out  of  Punta  Arenas,  Chile,  near  the  southern 
tip  of  South  America.  During  the  actual  experiment,  SAR 
imagery  with  coincident  aircraft  measurements  were  ac¬ 
quired  over  a  five-day  period.  By  applying  the  algorithm 
described  in  this  paper,  we  were  able  to  generate  esti¬ 
mates  of  two-dimensional  ocean  height-variance  spectra 
that  closely  agree  with  measurements  made  by  both  the 
ROWS  and  SCR  [12]. 

The  general  strategy  used  in  this  work  has  been  to  de¬ 
velop  algorithms  to  convert  SAR  imagery  into  height-var¬ 
iance  spectra  based  on  both  analytic  models  and  experi¬ 
ence  in  dealing  with  speckle  noise  inherent  in  SAR 
imagery.  The  effectiveness  of  the  algorithms  are  then  ver¬ 
ified  by  the  numerical  simulations  and  comparisons  with 
independently  measured  wave  spectra. 

II.  SAR  Wave  Imaging  Mechanisms 


effects  are  environmentally  dependent  and  generally  not 
well  understood  [  14]— [  16] . 

The  term  “radar  modulation  transfer  function”  is  used 
to  define  the  relationship  between  ocean  surface  wave 
height  and  the  radar  return  [7] .  The  effect  of  local  slope 
variations  is  described  by  the  tilt  modulation  transfer 
function 


R,  =  ikr  — 


(1) 


where  kr  is  the  range  component  of  ocean  wavenumber, 
<j0  is  the  radar  cross  section  of  the  surface,  and  9  is  the 
incidence  angle  [7].  Note  that  the  phase  of  this  transfer 
function  is  ±90°  with  respect  to  the  long  wave  crest,  de¬ 
pending  on  the  propagation  direction  of  the  wave  (i.e  , 
the  sign  of  kr).  The  magnitude  of  R,  can  be  estimated  by 
using  the  expression  for  <jo(0)  obtained  from  the  Bragg 
scattering  model  [13],  i.e., 

<to(0)  =  16ir*;5g2(0)  S(2*o  sin  9)  (2) 

where  *o  is  the  radar  wavenumber,  S(k)  is  the  wave  height- 
variance  spectral  density,  which  is  evaluated  at  the  Bragg 
wavenumber,  and  g(9)  =  I  ±  sin2  9  for  a  perfectly  con¬ 
ducting  surface  where  the  upper  sign  ( + )  holds  for  ver¬ 
tical  polarization  and  the  lower  (-)  sign  for  horizontal 
polarization.  Assuming  a  Phillips  spectrum  for  the  short 
waves,  i.-  .,  S(k)  =  Ak~ 4,  where  k  is  surface  wavenumber 
magnitude,  the  tilt  modulation  transfer  function  may  be 
approximated  by 


R,  =  ikr 


4  cot  9 
1  ±  sin2  9 


At  typical  incidence  angles  (20°  to  60°),  a  SAR  interacts 
with  the  ocean  surface  through  a  “Bragg”  resonance 
mechanism  with  short  ocean  surface  waves  [13].  The  ra¬ 
dar  cross  section  of  the  surface  is  dependent  on  the  local 
incidence  angle  and  the  height-variance  spectral  density 
of  the  surface  at  the  Bragg  wavenumber.  In  the  case  of  L- 
band  SAR,  like  Seasat  and  SIR-B,  the  Bragg  waves  are  in 
the  decimeter  wavelength  regime.  The  modulation  of  ra¬ 
dar  cross  section  by  long  ( >  50  m)  waves  renders  the  long 
waves  visible  in  real  aperture  radar  as  well  as  SAR  im¬ 
agery.  In  addition.  SAR  imagery  is  influenced  by  the  mo¬ 
tion  of  the  surface,  with  the  result  that  long  waves  may 
be  either  more  or  less  visible  in  a  SAR  image  than  in  a 
real  aperture  radar  image  with  the  same  nominal  resolu¬ 
tion. 

The  modulation  of  the  backscattered  power  that  is 
caused  by  changes  in  the  local  incidence  angle  is  desig¬ 
nated  as  the  tilt  effect.  The  modulation  in  power  caused 
by  changes  in  the  spectral  density  at  the  Bragg  wavenum¬ 
ber  across  the  profile  of  the  long  waves  is  called  the  hy¬ 
drodynamic  effect.  Changes  in  Bragg  wave  spectral  den¬ 
sity  are  due  to  interactions  of  these  short  surface  waves 
with  the  local  currents  induced  by  the  longer  waves,  and 
also  may  include  other  effects  such  as  those  due  to  vari¬ 
ations  of  wind  speed  over  the  long  wave  profile.  These 


over  the  range  of  incidence  angles  for  which  Bragg  scat¬ 
tering  is  valid. 

A  similar  transfer  function  can  be  defined  for  hydro- 
dynamic  modulation.  Theoretical  calculations  of  this 
transfer  function,  based  on  interactions  of  the  Bragg  waves 
with  the  orbital  currents  associated  with  longer  waves,  in¬ 
dicate  a  magnitude  of  4.5*2/*  [7],  assuming  a  Phillips 
spectrum  for  the  short  waves  and  the  relaxation  constant 
is  zero.  This  is  roughly  a  factor  of  two  smaller  than  the 
magnitude  of  the  tilt  modulation  transfer  function  at  9  = 
25°. 

The  effects  of  surface  motions  on  the  SAR  imaging  pro¬ 
cess  may  be  divided  into  two  categories.  First,  the  orbital 
velocities  associated  with  a  given  wave  component  cause 
periodic  displacements  within  the  SAR  image  that  result 
in  a  periodic  modulation  of  the  image  intensity.  If  these 
displacements  are  small  compared  to  the  azimuthal  com¬ 
ponent  of  wavelength,  the  image  modulations  caused  by 
this  “velocity  bunching"  effect  may  be  described  by  a 
transfer  function  whose  magnitude  is  given  by 

ft 

/?„  =  -  ^*0(stn2  9  sin2  ©  +  cos2  9)' 2  (4) 

where  R  is  the  range  between  the  SAR  platform  and  the 
surface,  V  is  the  platform  ground  track  velocity,  w  = 
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S4J 


•Jgk  is  the  radian  frequency  of  the  ocean  wave,  ka  is  the 
azimuth  wavenumber,  and  sin  <t>  =  k/k  [7].  If  the  inci¬ 
dence  angle  9  is  small,  this  transfer  function  reduces  to 


Rv  =  —  -Jgk  ka  cos  9. 


(5) 


The  phase  of  the  velocity  bunching  modulation  transfer 
function,  to  this  level  of  approximation,  is  either  zero  or 
180°,  depending  on  the  sign  of  ka. 

At  near  nadir  incidence  ( ?30°),  the  magnitude  of  the 
tilt  modulation  transfer  function  dominates  the  hydrody¬ 
namic  modulation  transfer  function.  Therefore,  we  have 
chosen  to  approximate  SAR  wave  imaging  at  near  nadir 
incidence  by  the  addition  of  the  tilt  and  velocity  bunching 
modulation  terms.  The  ultimate  justification  of  the  utility 
of  this  approximation  rests  in  its  effectiveness  in  estimat¬ 
ing  slope-  and  height-variance  spectra  from  SAR  im¬ 
agery. 

The  effects  of  velocity  bunching  and  tilt  modulation 
may  be  combined  by  adding  these  complex  modulation 
transfer  functions.  Since  these  transfer  functions  are  or¬ 
thogonal  (i.e.,  their  phases  differ  by  ±90°),  the  squared 
magnitude  of  the  combined  transfer  function  is  the  sum 
of  the  square  magnitude  of  each,  i.e.. 


where 


and 


i .  *>  ■* .  ■» 

flsAR  =  H'ka  +  m'k; 


m 


4  cot  9 
1  ±  sin2  9 


(6) 

(7) 


(8) 


When  the  wave  amplitude  becomes  sufficiently  large, 
so  that  the  maximum  scatterer  displacements  caused  by 
the  orbital  motion  are  larger  than  roughly  one-quarter  of 
the  azimuth  wavelength,  the  velocity  bunching  process 
becomes  nonlinear  in  the  sense  that  a  further  increase  in 
wave  amplitude  or  wave  height-variance  spectral  density 
does  not  result  in  a  proportional  increase  in  image  inten¬ 
sity-variance  spectral  density  at  the  same  wavenumber. 

This  effect,  along  with  the  effects  of  the  variations  in 
the  surface  velocity  within  each  resolution  cell  and  during 
the  SAR  integration  time,  form  a  second  category  of  mo¬ 
tion  effects.  These  effects  cause  a  departure  from  the  be¬ 
havior  predicted  by  (6).  Numerical  simulations  [9]  indi¬ 
cate  that  these  effects  may  be  accounted  for  by  multiplying 
the  transfer  function  given  in  (6)  by  an  azimuthal  weight¬ 
ing  or  falloff  function.  In  effect,  the  image  displacements, 
due  to  ocean  surface  motions  that  are  uncorrelated  with 
the  dominant  wave  system,  smear  the  image  in  azimuth. 
The  azimuth  falloff  function,  therefore,  can  be  thought  of 
as  a  low-pass  filter.  In  Section  V,  this  will  be  discussed 
in  more  detail. 
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Fig.  1.  Outline  of  procedure  to  estimate  ocean  wave  slope-  and  height- 
variance  spectra  from  SAR  imagery. 


III.  The  Estimation  of  Wave  Height  Spectra 
from  SAR  Imagery 

The  procedure  used  to  estimate  two-dimensional  slope- 
and  height-variance  spectra  from  SAR  imagery  can  be 
broken  into  a  five-level  process  outlined  in  Fig.  1.  It  is 
important  to  begin  with  SAR  images  that  are  both  radi- 
ometrically  corrected  and  geometrically  rectified.  For  the 
SIR-B  mission,  such  imagery  has  been  provided  to  us  by 
the  Jet  Propulsion  Laboratory. 

From  SAR  intensity  imagery,  a  single  SAR  image  frame 
can  be  extracted.  Such  a  frame  is  comprised  of  5 12  x  512 
image  pixels.  Since  each  pixel  represents  a  12.5  m  x  12.5 
m  area,  the  entire  image  frame  corresponds  to  a  6.4  km 
x  6.4  km  patch  on  the  ocean  surface.  This  frame  size 
provides  a  sufficiently  large  area  that  at  least  10  cycles  of 
very  long  surface  waves,  up  to  640  m  in  length,  can  be 
included  in  a  single  image  frame.  They  are  also  small 
enough  that  the  ocean  can  be  reasonably  assumed  homo¬ 
geneous  within  a  frame  [12], 

The  mean  image  intensity  is  subtracted  from  the  image 
frame  and  divided  into  the  result.  This  image  of  fractional 
modulation  is  then  Fourier  transformed  and  squared  to 
produce  a  level  1  image  intensity-variance  spectrum  as  a 
function  of  azimuth  and  range  wavenumber  kr). 

Because  a  SAR  image  has  finite  resolution  in  both  the 
azimuth  and  range  directions,  spectral  response  in  a  level 
1  spectrum  falls  off  at  high  azimuth  and  range  wavenum¬ 
bers.  The  estimation  of  and  correction  for  this  falloff  in 
response  has  been  reported  elsewhere  [2],  [17],  We  have 
designated  this  correction  as  the  stationary  response  cor¬ 
rection  to  distinguish  this  spectral  resoonse  function  from 
the  falloff  in  spectral  response  caused  b,  ocean  surface 
motion,  which  will  be  discussed  later.  The  stationary  re- 
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sponse  function  is  estimated  by  examining  SAR  imagery 
of  homogeneous,  spectrally  white  scenes.  The  falloff  is  lit 
to  a  model  function  which  is  then  inverted  and  used  as  a 
multiplicative  correction  to  the  image  spectrum.  After  this 
correction,  the  result  is  a  level  2  spectrum  S2(ka,  kr). 

Because  this  level  2  spectrum  has  only  2  degrees  of 
freedom,  the  value  of  the  spectrum  at  each  two-dimen¬ 
sional  wavenumber  bin  is  a  very  noisy  representation  of 
the  underlying  vpectrum.  The  underlying  spectrum  is  the 
one  that  in  principle  would  be  the  limit  obtained  by  the 
ensemble  average  of  many  spectral  measurements.  Spe- 
eitican « .  for  only  2  degrees  of  freedom,  the  rms  fractional 
error  in  each  ->pectral  estimate  is  100  percent. 

To  alleviate  this  problem,  the  level  2  spectrum  is  con¬ 
volved  with  a  Gaussian-shaped  kernel  sufficiently  wide 
that  each  spectral  estimate  in  the  level  3  spectrum  has 
about  300  degrees  of  freedom  w  ith  an  associated  rms  error 
of  only  8.5  percent  (2].  Of  course,  in  the  level  3  spectrum 
adjacent  spectral  estimates  are  no  longer  statistically  in¬ 
dependent  causing  a  reduction  in  resolution  in  wavenum¬ 
ber  space. 

It  is  possible  at  this  point  to  further  reduce  the  uncer- 
taintv  by  averaging  level  3  spectra  from  adjacent  image 
frames.  Assuming  the  ocean  surface  is  spatially  homo¬ 
geneous  over  these  frames,  the  averaged  spectra  still  will 
have  greater  statistical  reliability.  The  decision  to  average 
adjacent  spectra  is  dependent  on  the  particular  oceano¬ 
graphic  situation  considered.  In  regions  of  the  ocean  far 
from  current  boundaries  and  high  winds,  it  might  be  rea¬ 
sonable  to  average  spectra  over  large  areas.  In  regions 
where  spectra  might  be  more  rapidly  spatially  evolving, 
further  averaging  spectra  would  not  be  justified. 

SAR  imagery  is  contaminated  by  the  presence  of 
speckle  noise.  As  the  number  of  incoherent  looks  aver¬ 
aged  to  form  a  SAR  image  increases,  the  level  of  noise  is 
reduced.  This  muliplicative  noise  also  will  contaminate 
spectra  of  the  SAR  imagery.  The  amount  of  noise  present 
in  the  spectrum  is  dependent  on  both  the  mean  image  in¬ 
tensity  and  the  image  intensity  variance  caused  by  waves 
in  the  image  [18]. 

In  addition  to  the  speckle  noise  present  in  both  the  SAR 
imagery  and  resultant  spectra,  there  appears  to  be  another 
source  of  broad-band  spectral  noise.  One  possible  source 
is  a  "spikiness”  in  the  radar  return,  which  would  produce 
a  strongly  peaked  spatial  autocorrelation  function  of  the 
radar  cross  section,  or  equivalently  a  broad  spatial  spec¬ 
trum  of  the  radar  cross  section.  Unlike  the  effects  of  spec¬ 
kle,  the  contribution  of  such  radar  signal  spikes  to  the 
SAR  image  spectrum  would  be  modified  by  the  motion- 
induced  azimuth  falloff  function.  Clearly,  more  research 
is  needed  to  understand  the  source  of  these  effects.  How¬ 
ever,  in  most  cases,  an  empirical  procedure  can  be  fol¬ 
lowed  to  separate  the  noise  floor  from  the  “true”  spec¬ 
trum. 

To  remove  the  effects  of  broad-band  noise  in  the  spec¬ 
tra.  we  generate  level  4  spectra  by  subtracting  a  constant 
from  all  spectral  values  in  the  level  3  spectrum.  Any  neg¬ 
ative  values  are  set  equal  to  zero.  The  constant  level  to 


be  subtracted  is  calculated  using  the  results  from  Goldfin- 
ger [ 1 8  j . 

As  mentioned  in  the  previous  section,  the  SAR  modu¬ 
lation  transfer  function  defines  the  relationship  between 
the  SAR  image  spectrum  and  the  ocean  height-variance 
spectrum.  Employing  (6)  we  relate  a  level  4  spectrum  to 
a  level  5  height-variance  spectrum  by 


S?(ka,  1^) 


Sj(ka,  kr)  _  Sj{ka,  kr) 
flsAR  m'k';  4-  fji:k2a 


(9) 


Ss(ka,  kr )  is  the  level  5  estimate  of  the  height-var¬ 
iance  spectrum  of  the  ocean  surface.  Similar,  the  slope- 
variance  spectrum  Sf(ka,  kr)  =  k:S<(ka,  kr )  can  be  esti¬ 
mated 


kr) 


St(ka,  kr)  k1  _  Sx(ka,  kr)  k: 
R2sar  m:k;  ~  n:k:a 


(10) 


Although  the  integral  over  wavenumber  space  of  level 
5  spectra  given  in  (9)  and  (10)  ought  to  yield  estimates  of 
the  total  height  or  slope  variance  of  the  ocean  surface, 
the  normalization  of  the  image  spectrum  is  not  yet  well 
understood.  As  a  more  reliable  procedure,  level  5  spectra 
can  be  normalized  to  height  variance  or  slope  variance 
measured  independently.  Eventually,  we  anticipate  that 
further  analysis  of  SIR-B  and  independently  measured 
wave  spectra  will  result  in  more  accurate  estimates  of  the 
total  modulation  transfer  function. 


TV.  Comparison  between  SAR  and  Independently 
Measured  Wave  Spectra 

The  instruments  aboard  the  NASA  P-3  aircraft,  which 
measured  wave  spectra  in  the  vicinity  of  SAR  imagery, 
allow  us  to  verify  the  above  procedure  for  estimating  wave 
spectra  from  SAR  imagery.  Fig.  2  shows  a  set  of  four 
spectra.  The  spectrum  on  the  upper  left  is  a  level  5  SAR 
slope  spectrum  from  SIR-B  imagery  off  the  coast  of  Chile 
on  Shuttle  orbit  91 .  Note  that  azimuth  wavenumber  direc¬ 
tion  is  parallel  to  and  the  range  wavenumber  direction  is 
perpendicular  to  the  spacecraft  flight  direction.  The  spec¬ 
trum  on  the  upper  right  is  a  Radar  Ocean  Wave  Spectrom  ¬ 
eter  (ROWS)  [11]  estimated  slope-variance  spectrum 
measured  from  the  P-3  aircraft  approximately  3  h  aftei  the 
shuttle  overpass  at  the  same  spot  in  the  ocean.  Note  the 
extreme  similarity  in  the  spectra.  Both  show  a  bimodal 
wave  system  with  peaks  in  each  spectrum  showing  similar 
relative  magnitude. 

The  peaks  in  the  SAR  spectra  from  orbit  91.  however, 
are  predominantly  range  traveling.  Azimuth  traveling 
waves  are  somewhat  more  problematic  for  measurement 
by  a  SAR.  The  spectrum  in  the  lower  left  of  Fig.  2  is  a 
level  5  slope  spectrum  from  SIR-B  orbit  106  off  the  coast 
of  Chile.  The  peak  energy  system  is  almost  exactly  azi¬ 
muth  traveling  with  a  wavelength  of  approximately  350 
m.  A  secondary  wave  system  also  is  propagating  in  the 
range  direction.  The  lower  right-hand  comer  of  Fig.  2 
shows  the  corresponding  slope-variance  spectrum  from  the 
ROWS  Although  the  spectra  are  separated  by  about  50 
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Fig.  2.  Companson  of  ROWS  and  SAR  estimates  of  slop«-vanance  spec¬ 
tra.  (a)  Level  5  SAR  slope  spectrum  for  orbit  91.  time:  0230  October 
1 1 .  1984  (GMT),  location:  55.5°S.  82.5°W  (b)  ROWS  slope  spectrum, 
time:  0330  October  II.  1984  (GMT),  location:  55.3°S.  8l.2°W.  (c) 
Level  5  SAR  slope  spectrum  for  orbit  106.  time:  0045  October  12.  1984 
(GMT),  location:  56.00°S.  82. 0°W.  (d)  ROWS  slope  spectrum,  time: 
0400  October  12.  1984  (GMT),  location:  55.5°S.  81.5°W.  Extracted 
from  ( 12]. 

km,  they  both  confirm  the  existence  of  both  wave  sys¬ 
tems. 

The  comparison  of  level  5  SAR  slope  spectra  with  ac¬ 
tual  measurements  of  slope-variance  spectra  from  the 
P-3  aircraft  strongly  support  the  conclusion  that  level  5 
SAR  slope  spectra  are  reasonable  estimates  of  the  ocean 
slope-variance  spectrum. 

Fig.  3  shows  a  comparison  of  the  level  5  height  spectra 
from  SIR-B  with  corresponding  height-variance  spectra 
measured  by  the  Surface  Contour  Radar  (SCR)  [10]  for 
orbit  106.  The  SCR  was  not  operating  during  orbit  91. 
Again  note  the  similarity  between  the  spectra.  Of  course, 
the  k 2  conversion  between  the  slope  and  height  spectra 
tends  to  suppress  high-frequency  peaks  apparent  in  the 
slope  spectra.  For  this  data  set,  the  level  5  SAR  height 
spectra  are  highly  correlated  to  independently  measured 
wave  height  spectra. 

In  our  discussion  thus  far,  we  have  not  considered  any 
correction  for  the  fallofF  in  azimuth  response  associated 
with  ocean  surface  motion  [6],  (19).  A  reliable  systematic 
procedure  to  perform  this  correction  is  still  under  devel¬ 
opment.  As  was  discussed  in  Section  II  and  will  be  pointed 
out  in  Section  V,  the  effect  of  ocean  surface  motion  on  a 
SAR  image  is  proportional  to  the  R/V  ratio  of  the  SAR 
platform.  Because  SIR-B  orbited  at  a  230-km  altitude  as 
opposed  to  the  800-km  altitude  of  Seasat,  the  R/V  ratio 
for  SIR-B  is  approximately  35  s  as  compared  to  the  Seasat 
R/V  ratio  of  120  s.  The  next  section  will  demonstrate  that 
for  all  but  the  most  extreme  sea  states,  the  azimuth  falloff 


Fig.  3.  Companson  of  SCR  and  SAR  estimates  of  height-vanance  spectra 
(a)  Level  5  SAR  height  spectrum  for  orbit  106,  time:  045  October  12. 
1984  (GMT),  location:  56.0°S,  82.0°W.  (b)  SCR  height  spectrum  ob¬ 
tained  from  composite  of  three  flight  headings  (209°.  256°.  and  304° 
with  respect  to  North)  on  October  12,  1984.  Extracted  from  (12). 

only  minimally  affects  SIR-B  spectra  out  to  a  wavenum¬ 
ber  of  27r/(  100  m)  or  0.0628  rad/m.  which  are  the  outer 
limits  of  the  spectra  displayed  in  Figs.  2  and  3.  Therefore, 
we  leave  this  correction  for  the  future. 


V.  The  Limitations  of  SAR  Imagery  in  High  Sea 
States 

A  SAR  relies  on  Doppler  information  to  achieve  its  high 
azimuth  resolution.  Ocean  surface  motion,  which  alters 
this  Doppler  information,  serves  to  complicate  the  imag¬ 
ing  process.  In  Section  II,  the  advection  of  surface  scat¬ 
tered  by  the  orbital  velocity  of  ocean  surface  waves  is 
shown  to  be  the  source  of  the  velocity  bunching  wave  im¬ 
aging  mechanism  for  azimuth  traveling  waves. 

Velocity  bunching  serves  to  make  azimuth  traveling 
waves  visible  in  SAR  imagery  because  surface  velocities 
are  modulated  on  the  spatial  scale  of  the  long  ocean  wave 
and  are  highly  correlated  with  the  ocean  surface  slope  of 
that  wave.  However,  the  orbital  velocities  imparted  to 
surface  scattered  are  not  all  spatially  correlated.  The  mo¬ 
tion  of  scattered  induced  by  waves  of  wavelengths  smaller 
than  or  of  the  same  order  of  the  SAR  spatial  resolution 
certainly  do  not  contribute  to  the  imaging  of  waves.  Even 
at  large  spatial  scales,  some  fraction  of  the  velocity  vari¬ 
ance  on  the  ocean  surface  is  sufficiently  random  that  it  is 
destructive  to  the  image  forming  process.  Depending  on 
the  particular  SAR  configuration  and  the  sea  state.  SAR 
spectral  response  falls  rapidly  with  azimuth  wavenumber. 


B-7 


548 


IEEE  TRANSACTIONS  ON  GEOSCIENCE  AND  REMOTE  SENSING.  VOL.  GE-24.  NO  4.  JULY  I486 


At  sufficiently  large  azimuth  wavenumbers,  no  construc¬ 
tive  imaging  is  possible. 

In  this  section  we  will  endeavor  to  provide  a  heuristic 
explanation  for  this  lack  of  azimuth  response.  The  results 
of  numerical  simulations  and  measurements  made  from 
actual  SAR  spectra  are  presented  and  interpreted  in  the 
context  of  this  heuristic  explanation. 

Essentially  we  will  model  the  effect  of  ocean  surface 
motion  as  a  linear  low-pass  filter  applied  to  SAR  imagery. 
For  a  scatterer  with  a  velocity  vr  parallel  to  the  SAR  radar 
beam,  the  scatterer  will  be  displaced  in  azimuth  in  the 
SAR  image  by  an  amount  dx  given  by 


If  we  assume  that  the  spatially  random  portion  of  sur¬ 
face  velocity  variance  radial  to  the  SAR  is  described  by  a 
Gaussian  probability  distribution  with  rms  velocity  v0, 
then  the  azimuth  position  shifts  of  the  surface  scatterers 
also  will  be  Gaussian  distributed  with  an  rms  position  shift 
given  by 


These  random  position  shifts  can  be  thought  of  as 
smearing  the  SAR  image  in  the  azimuth  direction  or 
equivalently  of  convolving  the  unsmeared  SAR  image 
with  an  impulse  response  function  h{x)  given  by 

1  r  ,r:  1 

>K.c)  =  — 4—  exp  ?l.  (13) 

'J2vx,  l 

This  convolution  in  the  SAR  image  spectrum  is  equivalent 
to  multiplying  the  SAR  image  spectrum  of  the  unsmeared 
image  by 

r  I'-  1 

H(ka)  =  exp  1  1  (14) 

1  ■ 

where  k,  is  given  by 


V 


This  A:,  is  the  azimuth  wavenumber  where  spectral  re¬ 
sponse  falls  to  0.61  of  its  maximum  value.  We  define  X, 
as  the  wavelength  associated  with  Ac,  where  X,  is  given  by 

r-  R 

X,  =  2ir  v2  -  vT  (16) 

The  quantity  X,  is  not  literally  the  azimuth  "cutoff" 
wavelength,  since  azimuth  traveling  waves  somewhat 
shorter  than  X,can  be  imaged.  Experimental  observations 
and  numerical  simulations  of  SAR  imagery  discussed  be¬ 
low  indicate  that,  for  the  spectral  smoothing  described  in 
Section  III.  the  azimuth  cutoff  wavelength  is  about  one- 
half  X,. 

Attempts  to  estimate  u,  have  been  made  by  Beal  et  al. 
[1]  and  Thomas  [20) .  [21].  Clearly,  u,  is.  in  general,  de¬ 


pendent  upon  both  the  ocean  height-variance  spectrum  and 
relative  SAR  geometry.  Since  Seasat  had  a  nadir  inci¬ 
dence  angle  of  20°  and  wave  imagery  of  .ained  from 
SIR-B  is  generally  at  similar  nadir  angles,  we  can  reason¬ 
ably  neglect  geometric  considerations  and  consider  only 
the  vertical  component  of  orbital  velocity  as  being  parallel 
to  the  radar  beam.  At  20°  from  nadir,  results  should  be 
accurate  to  about  10  percent. 

If  S(uj)  is  defined  as  the  one-dimensional  height-van- 
ance  frequency  spectrum,  then  uj:S(w)  is  the  velocity-var¬ 
iance  frequency  spectrum.  For  deep  water  waves,  the  ver¬ 
tical  component  of  orbital  velocity  is  equivalent  to  the 
horizontal  component  and  is  given  by 

r  1 1  - 

V„  =  \  (lui  i  i  Pt 

__  J  Wl 

The  velocities  contributing  to  (17)  may  be  considered 
as  the  velocities  of  scattenng  patches  or  facets  in  a  two- 
scale  model  [22].  Thus,  the  upper  limits  of  integration  oj; 
should  correspond  to  the  cutoff  wavenumber  in  the  two- 
scale  model.  However,  the  contributions  from  waves 
shorter  than  this  cutoff  is  usually  very  small.  Therefore, 
for  the  purpose  of  these  calculations,  the  upper  limits  ot 
integration  ^  can  be  taken  to  be  infinity 

Tucker  [20)  has  suggested  that  only  surface  velocities 
induced  by  waves  with  wavelengths  less  than  twice  the 
SAR  resolution  limit  contribute  to  the  smearing  Selecting 
the  lower  limit  of  integration  to  correspond  to  the  fre¬ 
quency  of  a  wave  having  a  wavelength  X  ,  equal  to  twice 
the  resolution  of  the  SAR.  and  assuming  that  S imi  has  the 
general  form  of  the  Pierson-Moskowitz  spectrum  [23], 
the  appropriate  rms  velocity  that  contributes  to  smeanng 
is  given  by 

u„  =  (0.442  m1  "s>  H\  :[erf  (2.86  x  I0':  X,  : 

i  IS) 

where  Hf  is  significant  wave  height.  The  corresponding 
X,  is 

X;  =  (3.93  m' -/s)  *  H\  :[erf  (2.86  x  10‘:X, 

i  I9i 

If  we  assume  that  velocities  at  all  spatial  scales  contribute 
to  azimuth  smeanng,  i.e..  X,  —  oo  the  worst  case  as¬ 
sumption.  ( 19)  reduces  to 

X®  =  (3.93  m'  :'s)  p  H\ :  1 20) 

Regardless  of  what  fraction  of  the  total  surface  velocity 
vanance  contnbutes  to  azimuth  smeanng,  ( 1 4)  and  (20) 
make  clear  the  cntical  dependence  of  azimuth  smeanng 
on  the  R/V  ratio  of  the  SAR  platform. 

These  results  may  be  compared  with  numencal  simu¬ 
lations  of  SAR  image  spectra.  Fig.  4  shows  the  total 


B-8 


monaldo  and  lyzenga  wave  slope-  and  height-variance  spectra 


549 


Fig  4  Ratio  of  image  spectrum  to  height  spectrum  along  azimuth  wave¬ 
number  axis  for  S  T  »  30  s  obtained  from  numerical  simulation  of  SAR 
image  spectra  First  number  associated  with  each  curve  indicates  domi¬ 
nant  wavelength  i in  meters),  the  second  indicates  significant  wave  height 
•  in  meters)  Extracted  from  (9). 


SIGNIFICANT  WAVE  HEIGHT  (meters) 

Fig.  5.  A  companson  of  estimated  from  (19)  and  (20)  and  the  Seasat 
and  SIR-B  R/V  ratios.  Solid  lines  refer  to  (20),  dashed  lines  to  (19)  as¬ 
suming  =  50  m.  Crosses  marked  X„'s  estimated  from  actual  Seasat 
SAR  image  spectra.  Circles  marked  \„'s  estimated  from  SIR-B  SAR  im¬ 
age  spectra.  The  numbers  refer  to  orbit  numbers,  the  names  refer  to  hur¬ 
ricanes  in  the  vicinity  of  SAR  imagery. 


table  i 

Comparison  of  Predicted  Using*  Numerical  Simulation  and  (19) 
and  1 20)  for  Various  Significant  Wave  Heights  and 
Dominant  Wavelengths 
i  R/V  •  30  s.) 


Significant  , 
“eignt 
meters 

dominant 
««*«( engt* 
meters  1 

3  l 

Simulations 
.meters,  I 

i* 

3 

Eauatton  ’9 
(meters) 

|  ‘5 

|  Equation  20 
j  (meters) 

2 

*30 

’92 

•39 

!  '67 

:oo 

235 

'47 

!  ns 

i 

400 

105 

1 50 

;  334 

-:o 

*47 

68 

82 

; 

!3C 

'27 

71 

116 

4 

400 

58 

74 

’67 

1 

transfer  function  obtained  from  the  simulation  model  de¬ 
scribed  in  (9)  for  live  different  wave  conditions.  The  cases 
for  which  the  significant  wave  height  is  equal  to  two  per¬ 
cent  of  the  dominant  wavelength  correspond  to  fully  de¬ 
veloped  Pierson-Moskowitz  spectra.  The  remaining  cases 
use  the  same  spectral  shape,  but  have  a  total  significant 
wave  height  reduced  by  a  factor  of  two. 

The  transfer  functions  in  Fig.  4,  shown  for  the  azimuth 
wavenumber  direction,  may  be  fit  quite  closely  by  the 
product  of  ( 19)  and  (14)  with  the  values  of  k ,  shown  in 
Table  I.  The  corresponding  X7's  also  are  shown  in  Table 
I  along  with  X '  and  Xf  computed  from  (19)  and  (20).  re¬ 
spectively. 

The  results  for  the  fully  developed  spectra  shown  in  the 
Table  indicate  that  the  azimuth  falloff  functions  obtained 
from  the  numerical  simulations  are  closer  to  the  predic¬ 
tions  of  (20).  i.e. .  the  integration  of  velocity  variance  over 
the  entire  spectrum,  than  the  predictions  using  (19).  That 
is.  the  calculations  agree  reasonably  well  with  the  numer¬ 


ical  simulations  based  on  the  Pierson-Moskowitz  spec¬ 
trum. 

The  agreement  is  also  good  for  the  cases  where  the  sig¬ 
nificant  wave  height  of  the  Pierson-Moskowitz  spectrum 
is  reduced  by  a  factor  of  two,  except  for  the  shortest 
wavelength  (100  m)  case  considered.  In  this  case,  the  de¬ 
parture  from  a  “linear”  transfer  function  appears  to  be 
caused  by  noniinearities  in  the  imaging  of  the  dominant 
wave  itself,  whereas  in  the  other  cases,  the  smearing  ef¬ 
fects  from  random  surface  motion  are  most  important. 
That  is.  for  this  particular  case,  the  treatment  of  the  ef¬ 
fects  of  ocean  surface  motion  as  a  low-pass  filter  in  a  lin¬ 
ear  system  starts  to  break  down.  Except  for  this  case,  the 
SAR  imaging  simulation  supports  the  conclusion  that  the 
effects  of  motion  generally  can  be  treated  as  a  low-pass 
filter  on  a  linear  system. 

To  examine  how  well  the  analytic  and  numerical  mod¬ 
eling  describe  the  azimuth  falloff  problem,  we  estimated 
the  falloff  observed  in  a  series  of  SAR  spectra  from  both 
Seasat  and  SIR-B.  The  estimate  of  the  azimuth  falloff  was 
performed  by  excluding  from  consideration  the  regions  in 
each  spectrum  having  apparent  concentrations  of  wave 
spectral  energy.  The  remaining  regions  of  the  spectrum 
exhibited  a  gradual  spectral  energy  falloff  in  the  azimuth 
direction  resting  on  a  relatively  low  flat  platform  of  spec¬ 
tral  energy  density. 

After  subtraction  of  this  platform,  the  spectral  falloff 
was  fit,  in  a  least  squares  sense,  to  the  Gaussian-shaped 
azimuth  falloff  as  suggested  by  (14).  The  resulting  fit 
yielded  the  parameter  k„  or  alternatively  X„. 

Fig.  5  compares  estimates  of  azimuth  falloff  for  various 
SAR  image  spectra  from  Seasat  and  SIR-B  with  expected 
values  of  X„.  The  solid  lines  represent  the  limiting  case, 
where  all  velocities  contribute  to  azimuth  smearing, 
yielding  Xf  proportional  to  H J,:.  for  the  Seasat  and  SIR- 
B  R/V  ratios.  The  dashed  lines,  asymptotically  indepen¬ 
dent  of  significant  wave  height,  represent  the  predicted 
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X»  assuming  only  velocities  on  spatial  scales  less  than  50 
m  (twice  the  SAR  spatial  resolution  for  Seasat  and  SIR- 
B). 

Note  that,  in  general,  the  observed  azimuth  falloffs  are 
significantly  below  values  predicted  in  either  (19)  or  (20). 
Assuming  that  velocities  on  spatial  scales  above  50mu 
not  contribute  to  azimuth  smearing  suggests  a  maximum 
X„,  asymptotically  approached  at  high  sea  states.  The  X, 
at  high  sea  states,  particularly  for  SIR-B  spectra,  appear 
to  approach  a  constant  at  high  sea  states ,  but  the  data  are 
not  conclusive. 

These  numerical  and  analytic  model  estimates  assume 
a  broad-band  fully  developed  sea.  In  SAR  image  spectra, 
we  rarely  observed  fully  developed  seas.  Much  of  the  ve¬ 
locity  variance  is  thus  highly  correlated  to  the  swell  wave 
slope  and  therefore  does  not  contribute  to  azimuth  smear¬ 
ing.  Thus,  the  extent  of  azimuth  smearing  is  somewhat 
less  than  expected  purely  on  the  basis  of  significant  wave 
height  and  (19)  and  (20),  and  the  simulations.  The  general 
leveling  off  of  X  „  at  high  sea  states  in  the  data  may  simply 
indicate  that  the  ocean  is  rarely  fully  developed  at  these 
high  sea  states. 

Equations  (19)  and  (20),  and  the  simulation  results 
when  using  broad-band  input  spectra  therefore  can  be  in¬ 
terpreted  as  upper  bounds  to  X„.  In  all  cases  observed,  X„ 
has  been  found  to  be  significantly  lower  although  the  de¬ 
pendence  on  R/V  is  corrected.  For  the  case  of  SIR-B,  even 
in  the  immediate  region  of  Hurricane  Josephine  whose 
significant  wave  height  was  estimated  at  8  m,  X„  was  on 
the  order  of  100  m.  For  Seasat  SAR  image  spectra  from 
Hurricane  Iva,  having  an  estimated  significant  wave  height 
of  6.5  m,  the  estimated  from  the  SAR  image  spectrum 
was  much  larger,  about  300  m. 

VI.  Conclusions 

A  procedure  has  been  proposed  to  estimate  ocean  slope- 
and  height-variance  spectra  from  SAR  imagery.  Results, 
based  on  numerical  simulation  of  SAR  image  spectra  and 
comparisons  between  SAR  image  spectra  from  SIR-B  and 
spectra  measured  by  independent  means,  suggest  that  SAR 
image  spectra,  when  multiplied  by  the  appropriate  linear 
transfer  functions,  produces  accurate  estimates  of  both  the 
slope-  and  height- variance  spectra.  The  resulting  SAR 
spectra  closely  agree  with  independent  measures  of  either 
the  ocean  wave  slope  and  height- variance  spectra. 

The  limitations  of  SAR  image  spectra  at  high  azimuth 
wavenumbers  in  high  sea-state  conditions  is  directly  pro¬ 
portional  to  the  R/V  ratio  of  the  SAR  platform.  For  R/V 
ratios  of  35  s,  typical  of  shuttle  altitudes,  significant  azi¬ 
muth  spectral  response  extends  to  ocean  wavelengths  of 
at  least  100  m  in  the  highest  expected  sea  states. 
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